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PREFACE 


The  David  W.  Taylor  Lectures  were  conceived  to  honor  our  founder  in 
recognition  of  his  many  contributions  to  naval  architecture  and  naval 
hydrodynamics.  Admiral  Taylor  had  a  keen  devotion  to  both  experimental 
evaluations  and  theoretical  analysis  for  the  solution  of  naval  problems. 

He  established  a  tradition  of  applied  scientific  research  at  the  "Model 
Basin"  which  has  been  carefully  nurtured  through  the  decades  and  which  we 
treasure  and  maintain  today.  It  is  in  this  spirit  that  we  have  invited 
Prof.  Allan  J.  Acosta  to  be  a  David  W.  Taylor  Lecturer. 

Dr.  Acosta  received  his  Bachelor  of  Science  Degree  from  the  California 
Institute  of  Technology  in  1945,  followed  by  a  Master  of  Science  Degree 
(1948)  and  the  Ph.D.  (Mechanical  Engineering)  in  1952.  For  the  following 
two  years,  he  was  Section  Chief  of  the  Hydrodynamics  Laboratory.  During 
the  period  1954-58,  Dr.  Acosta  was  Assistant  Professor  in  the  Mechanical 
Engineering  Department  at  the  California  Institute  of  Technology,  was 
promoted  to  Associate  Professor  in  1958,  and  since  1966  has  been  full 
professor.  Currently,  he  teaches  thermodynamics,  fluid  mechanics,  heat 
transfer,  and  advanced  turbomechanics.  His  research  areas  include  hydro¬ 
dynamics  of  cavitating  flows,  hydrofoils,  underwater  bodies,  turbomachines 
and  heat  transfer. 

Prof.  Acosta  lives  near  the  beach  and  is  an  avid  sailor.  He  is 
currently  completing  construction  of  a  sailboat  for  an  around-the-world 
voyage . 
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FOREWORD 


It  has  been  a  great  honor  to  deliver  the  fourth  in  this  series  of 
lectures  because  of  the  prestige  of  this  laboratory  and  the  distinguished 
lecturers  before  me.  The  previous  lectures  have  all  been  of  a  rather 
theoretical  nature;  1  have  been  asked  instead  to  emphasize  problems  from  a 
more  applied  point  of  view.  On  the  one  hand,  this  appears  to  be  an  easy 
course  to  follow,  since  many  important  areas  of  naval  hydrodynamics  lack 
the  precise  physical  description  required  for  detailed  mathematical 
analysis.  On  the  other  hand,  there  is  a  real  problem  of  selection,  since 
so  much  of  the  present  work  in  hydrodynamics  has  originated  here  at  the 
Naval  Ship  Research  and  Development  Center  (NSRDC) .  With  this  caveat  then, 
it  will  be  understandable  that  I  have  chosen  for  discussion  a  subject  of 
long  standing  interest  in  naval  hydrodynamics,  generally,  and  one  of  partic¬ 
ular  interest  to  this  laboratory,  namely,  cavitation.  Even  this  is  too 
broad  a  topic  for  only  a  few  lectures  and  I  have  restricted  myself  to  two 
topics;  these  are,  cavitation  inception  on  smooth  bodies  and  some  effects 
of  developed  cavitation  in  internal  flows.  There  are  many  instances  of 
cavitation  in  internal  flows.  The  testing  of  cavitating  bodies  in  water 
tunnels  is  one,  and  that  of  cavitation  in  a  ducted  propulsion  system  is 
another.  In  both  situations  the  flow  is  internally  confined  and  is,  there¬ 
fore,  geometrically  complex.  System  interactions  can  then  occur  which, 
particularly  for  unsteady  flows,  leads  to  great  difficulties  of  experi¬ 
mentation  and  interpretation. 

I  would  like  to  acknowledge  the  help  and  stimulating  discussions  with 
many  of  the  laboratory  staff.  These  are  too  many  to  mention  here,  but  I 
would  like  to  reserve  special  thanks  to  Justin  McCarthy,  William  Morgan, 
and  Mrs.  Shirley  Childers  for  making  my  lectureship  a  memorable  occasion. 
Finally,  I  would  like  to  express  my  gratitude  to  my  friend  and  former 
teacher,  Milton  Plesset,  for  his  witty  and  insightful  counsel. 
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CAVITATION  INCEPTION  AND  INTERNAL  FLOWS 
WITH  CAVITATION 

INTRODUCTION 

Cavitation  has  been  a  topic  of  unique  importance  in  naval  hydro¬ 
dynamics  since  the  introduction  of  high  speed  propellers.  The  interesting 
historical  survey  by  Johnson  (1972)*  shows  that  the  principal  features  of 
cavitation  as  a  source  of  noise,  vibration,  damage,  and  reduced  performance 
were  experienced  at  the  outset  in  ship  propulsion  in  a  pattern  that  has 
continued  to  be  repeated  in  the  most  current  applications.  These  un¬ 
desirable  features  of  cavitation  were  highlighted  in  Knapp's  well-known 
Clayton  Lecture  before  the  Institute  of  Mechanical  Engineers  (Knapp  1952) 
and  is  echoed  in  the  summary  volume  on  cavitation  (Knapp  et  al.  1970). 

The  inception  of  cavitation  announces  the  start  of  the  phase  change  that 
may  lead,  in  turn,  to  these  unwanted  effects  and  is  therefore  a  boundary 
between  fully  wetted  and  cavitating  modes  of  flow  that  is  very  desirable 
to  know.  But  the  onset  of  cavitation,  when  it  does  occur,  does  so  in  a 
variety  of  forms  which  may  differ  from  facility  to  facility  on  similar 
bodies.  Even  on  one  test  body,  the  various  appearances  of  cavitation  may 
change  with  flow  speed  and  liquid  environment. 

Thus,  this  inception  process  does  not  appear  to  be  a  well-defined, 
unique  kind  of  physical  event  for  which  there  are  clear-cut  simple  govern¬ 
ing  relations.  Indeed,  as  in  all  fluid  mechanics  the  presence  of  cavita¬ 
tion  requires  appeal  to  all  three  of  the  basic  conservation  laws  (mass, 

/ 

momentum,  energy) .  Cavitation  inception  appears  to  be  one  of  those  physical 
problems  which  does  not  yet  seem  well  posed  for  the  laboratory;  experi¬ 
mental  results  on  supposedly  similar  tests  show  wide  variability  and,  even 
for  one  specific  test  condition,  a  wide  "scatter"  is  often  observed.  This 
state  of  affairs  was  not  helped  much,  in  early  years,  by  the  absence  of  a 
detailed  scientific  description  of  the  flow  (with  some  exceptions)  at  the 
moment  of  cavitation  inception.  It  was  then  not  surprising  that  many  of 
the  theoretical  attempts  to  explain  the  "scaling  laws"  fell  short. 


*A  bibliographic  listing  of  references  is  given  on  page  93. 
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Cavitation  inception  as  a  separate  topic  of  study  has  seemed  then  to 
be  a  somewhat  "messy"  subject.  Yet,  it  is  an  important  one  in  fluid 
engineering,  and  recently  there  has  been  some  substantial  progress  by  the 
combined  efforts  of  workers  in  several  countries  including  this  Center, 
both  in  understanding  the  basic  flows  involved,  and  in  laboratory  tech¬ 
niques.  In  the  remaining  sections,  we  will  review,  briefly,  some  of  these 
findings  and  recent  experimental  work.  To  appreciate  the  point  of  view  of 
many  of  the  workers  in  cavitation  we  need  to  survey  the  highlights  of 
bubble  mechanics,  the  most  common  denominator  in  cavitation-inception 
scaling  theories.  At  the  root  of  all  phase  change  processes  in  cavitation 
is  a  "nucleating  source" — or  simply  a  cavitation  nucleus — since  technical 
fluids  do  not  exhibit  the  large  tensions  required  of  pure  liquids  prior  to 
cavitation.  These  nuclei,  long  the  subject  of  speculation,  are  now  be¬ 
coming  accessible  to  measurement,  although  there  is  still  an  observational 
gap  connecting  these  sources  of  cavitation  and  the  macroscopic  cavitation 
ultimately  seen  on  bodies.  Some  of  these  findings  are  reviewed  in  light 
of  current  laboratory  practice.  Finally,  we  discuss,  in  a  preliminary  way, 
some  of  the  special  problems  nosed  by  developed  cavitation  when  this  occurs 
in  confined  flows. 

CAVITATION  INCEPTION 
Preliminaries 

It  is  common  to  express  the  ambient  pressure  at  which  cavitation  occurs 
in  a  coefficient  form  called  the  "cavitation  number" 

p  -p  (T  ) 
roo  r  v  00 

o  = - - -  (1) 

2 

where  q  is  the  freestream  dynamic  pressure,  1/2  pU  ,  p  is  the  ambient 

o  00 

static  pressure  and  Pv(Too)  is  the  vapor  pressure  of  the  pure  liquid  at 
ambient  (i.e.,  freestream)  temperature.  At  the  same  time  the  static  pres¬ 
sure,  p,  on  a  body  is  put  into  dimensionless  form 
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(2) 


c 
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P“P„ 
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so  that  the  pressure  is 


p  =  pv(tJ  +q(o+cp) 


(3) 


When  the  fluid  pressure,  p,  is  less  than  the  vapor  pressure,  p^,  we  expect 
vapor  to  form  and,  hence,  we  may  expect  cavitation  to  begin.  Thus,  as  a 
first  guess,  we  estimate  that 


a  <  -  c 

^min 


(4) 


is  the  condition  for  the  onset  of  cavitation.  We  follow  Holl  (1969)  in 
calling  the  cavitation  that  occurs  when  Pm^n<Pv^T0O)  vaporous  cavitation  as 
we  expect  the  cavitation  process  to  be  one  of  the  formation  and  collapse 
of  cavitation  voids  filled  primarily  with  the  vapor  of  the  liquid. 
Observations  show  (Holl  ibid)  that  it  is  also  possible  for  Pmin  to  exceed 
the  vapor  pressure  at  inception.  This  can  happen  wherever  the  test  liquid 
contains  dissolved  gas  and  the  local  pressure  is  below  the  gas  saturation 
pressure  or  bubble  point.  Cavitation  inception,  when  so  dominated  by 
diffusion  of  the  gas  into  a  cavity  void,  is  termed  gaseous  cavitation  by 
Holl.  As  he  has  shown,  this  may  be  quite  common  in  hydrodynamic  appli¬ 
cations  (see  also  the  reviews  by  Acosta  and  Parkin  1975,  Morgan  and 
Peterson  1977).  A  further  type  of  cavitation  inception  termed  by  Holl 
"pseudo  cavitation,"  is  said  to  occur  when  entrained  bubbles  of  gas  or 
vapor  respond  quasi-statically  to  an  imposed  pressure  history.  This  fea¬ 
ture  of  cavitation  is  more  characteristic  of  bubbly  two-phase  flow  than  it 
is  of  cavitation  inception  and  will  not  be  considered  further  herein. 
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Scale  Effect — Early  Results 

Departures  from  the  rule 
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Pmin 
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for  the  cavitation  number  at  inception  are  called  "scale  effects."  Indeed, 
almost  from  the  very  start  of  experimentation  in  cavitation,  discrepancies 
from  this  rule  were  so  prevalent  that  this  became  an  important  research 
topic,  particularly  for  laboratory  work  where  the  discrepancies  were  more 
evident.  The  thrust  of  the  earlier  studies  was,  as  it  is  now,  to  obtain  a 
more  detailed  microscopic  picture  of  the  formation  of  cavitation  in  the 
hope  that  this  would  lead  to  a  better  physical  model  for  describing  the 
process. 

Among  the  very  most  influential  of  these  works  were  the  beautiful 
motion  pictures  of  growing  and  collapsing  bubbles  travelling  on  the  sur¬ 
face  of  a  body  of  revolution  made  by  R.T.  Knapp  of  Caltech  in  about  1946 
These  were  sensational  pictures  at  the  time  which,  when  coupled  with  the 
brilliantly  successful  analysis  of  Plesset  (see  Knapp  and  Hollander  1948, 
Plesset  1949),  firmly  established  bubble  mechanics  as  the  essential 
scientific  component  of  this  type  of  cavitation.  In  fact,  because  of  the 
time-delay  required  for  the  growth  of  a  travelling  bubble  to  visible  size, 
it  appeared  reasonable  to  expect  the  scale  effect  to  follow  from  this 
mechanism.  A  reproduction  of  one  of  Knapp's  photographs  of  these  travel¬ 
ling  bubbles  is  shown  in  Figure  1.  There,  then,  followed  in  the  early 
1950’ s  a  period  of  intense  activity,  both  experimental  and  theoretical,  de¬ 
voted  to  settling  the  cause  for  the  scale  effect.  To  a  latter  day  observer 
it  appears  that  bubble  mechanics,  already  of  vital  importance  in  under¬ 
water  acoustics,  became  the  paradigm  through  which  experiments  on  cavita¬ 
tion  inception  were  to  be  explained.  As  we  will  see,  growing  travelling 
bubbles  are  one  form  of  cavitation  characteristic  of  inception  but  they  are 
not  the  only  form,  nor  necessarily  the  dominant  form.  However,  at  this 
time  in  1950  these  differences  were  largely  overlooked  or  unknown. 
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Travelling  Bubble  Cavitation  on  a 
(Knapp  and  Hollander  1948) 


Figure  1 


A  good  account  of  cavitation  in  this  period  may  be  found  in  the  book 
"Cavitation"  (Knapp  et  al.  1970)  through  about  1969.  This,  together  with 
the  volume  "Cavitation  State  of  Knowledge,"  provides  an  excellent  summary 
to  that  time.  In  the  present  context,  some  of  the  events  of  this  period 
that  seem  significant  are  summarized  in  Table  1.  This  tabulation  stops 
in  1969  deliberately  to  emphasize  progress  made  since  then. 

TABLE  1  -  SOME  CHRONOLOGICAL  EVENTS  IN 
CAVITATION  INCEPTION 

1947  Knapp's  movies  of  growth  and  collapse  of  travelling 
bubbles 

1949  Plesset's  analysis  of  bubble  motion 

1953  Parkin  and  Kermeen  -  Boundary  Layer  Inception 

1955  Plesset  and  Zwick  -  Vapor  Bubble  Dynamics 

1956  Daily  and  Johnson  -  Shear  Flow  Inception 

1959  Ripkin  and  Killen  -  Nuclei  Measurements 

1962  van  der  Walle  -  "Stabilized"  Bubble  Cavitation 

1966  Johnson  and  Hsieh  -  Bubble  Screening 

1966  Schiebe  -  Cavitation  Occurrence  Counting 

1966  Lindgren  and  Johnsson  -  ITTC  Comparative  Tests 

1969  Holl  and  Kornhauser  -  Thermodynamic  Effects  on 
Scaling  Inception 

1969  "Cavitation  State  of  Knowledge" 


One  of  the  most  influential  events  of  this  tabulation  is  the  experi¬ 
mental  work  of  Parkin  and  Kermeen  (1953).  They  focus  more  sharply 
on  the  small  bubbles  seen  within  the  boundary  layer  of  smooth  bodies  as 
this  appeared  to  be  the  crucial  region  of  flow.  Their  impressive  photo¬ 
graphs  showed  that,  on  a  hemispher ically  nosed  body,  the  region  of  visible 
cavitation  was  preceded  by  a  region  of  "microscopic  cavitation  bubbles 
which  grow  in  the  boundary  layer."  Figure  2,  taken  from  their  report, 
shows  these  small  bubbles  upstream  of  the  more  readily  visible  macroscopic 
cavitation.  Also  shown  on  the  figure,  is  their  interpretation  of  the  flow 
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regime  on  the  body,  namely  an  attached  boundary  layer  in  which  these  small, 
essentially  stationary,  bubbles  grow,  they  conclude,  by  a  process  of  air 
diffusion.  Subsequently,  these  bubbles,  when  of  sufficient  size,  would  be 
"stripped  off"  to  "feed"  the  downstream  macroscopic  cavitation  seen  by  the 
eye.  It  will  be  noticed  that  at  the  location  of  the  observed  small  bubbles 
the  local  pressure  is  actually  above  the  fluid  vapor  pressure,  hence, 
vaporous  growth  would  be  impossible  there.  Presumably,  truly  microscopic 
bubbles  could  grow  further  upstream  near  the  minimum  pressure  point,  but 
these  were  not  observable  with  the  photographic  technique  used.  Parkin 
and  Kermeen  also  stressed  that  travelling-bubble  cavitation  did  not  occur 
in  their  tests  readily  because  the  new  facility  had  a  "resorber"  which 
prevented  comparatively  large  freestream  bubbles  (the  nuclei  for  the 
travelling  bubbles)  from  being  recirculated,  and  thus,  that  this  type  of 
cavitation  had  a  true  boundary-layer  origin. 

This  work  had  a  vital,  energizing  effect  on  cavitation  research  and 
the  physical  ideas  presented  were  used  as  the  basis  for  a  number  of 
theories  of  the  scaling  effect  (e.g.,  van  der  Walle  1962,  Holl  and 
Kornhauser  1969,  among  others).  These  theories  incorporated  the  idea  of  a 
microbubble  stabilized  on  the  surface  of  a  body  so  that  growth  by  gas 
diffusion  could  occur  until  some  critical  size  was  achieved  when  downstream 
transport  into  the  flow  would  take  place.  Further  refinements  in  bubble 
dynamics  accounting  for  thermodynamic  effects  (Plesset  and  Zwick  1955)  and 
growth  by  oscillating  pressure  fields,  i.e.,  "rectified  diffusion”  (Hsieh 
and  Plesset  1961),  did  not  provide  additional  mechanisms  sufficient  to 
explain  the  large  scale  effect.  Before  turning  to  some  of  these  data  we 
mention  two  additional  concepts  dealing  with  inception  on  smooth  surfaces. 

The  first  of  these,  due  to  Daily  and  Johnson,  pointed  out  that  the 
turbulent  fluctuations  in  fully-developed  pipe  flow  can  also  lead  to 
cavitation.  (This  observation  has  been  expanded  subsequently  by  Arndt  and 
others  to  deal  with  jet  cavitation,  but  this  is  not  the  subject  of  the 
present  work.)  In  the  meantime,  Ripken  and  Killen  were  able  to  demonstrate 
in  certain  types  of  cavitation  the  importance  of  the  undissolved  gas,  i.e., 
the  microbubble  "nuclei"  within  the  flow.  In  addition,  they  published  size 
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distributions  of  these  freestream  nuclei,  the  first,  to  our  knowledge,  to 
be  measured  in  a  flow  facility.  It  seemed  quite  plain  then,  that  these 
nuclei  were  the  ones  responsible  for  the  travelling  bubbles  photographed  by 
Knapp.  That  being  so,  it  would  appear  that  if  these  freestream  microbubbles 
were  diminished  in  number  and  size,  travelling-bubble  cavitation  would  be. 
similarly  diminished.  One  mechanism,  whereby  the  larger  freestream  bubbles 
could  be  selectively  removed  from  the  low  pressure  region  surrounding  a 
body,  was  proposed  by  Johnson  and  Hsieh;  in  their  concept  bubbles  migrate 
across  streamlines  under  the  action  of  the  pressure  field  around  the  body. 
Indeed,  they  show,  by  calculation,  that  the  population  of  larger  bubbles 
near  the  position  of  minimum  pressure  can  be  reduced,  and  that  from  this 
"screening"  effect  a  size-dependent  effect  on  travelling-bubble  growth 
with  a  given  freestream  nuclei  distribution  is  found.  Thus,  a  scale  effect 
is  produced  based  upon  the  concept  of  travelling  cavitation  bubbles 
originating  from  freestream  microbubbles.  At  about  the  same  time,  Schiebe 
(1966)  proposed  a  scheme  of  inception  measurement  based  upon  the  idea  of 
measuring  the  rate  of  these  travelling-bubble  cavitation  events,  a  concept 
which  has  subsequently  received  much  attention  and  development  at  the  St. 
Anthony  Falls  Hydraulic  Laboratory  and  elsewhere. 

This  early  development  of  cavitation-inception  observation,  develop¬ 
ment  of  basic  theory  and  phenomenological  tests,  may  be  said  to  have 
culminated  in  the  round  robin  tests  on  a  standard  body  sponsored  by  the 
International  Towing  Tank  Conference  in  1966.  But  before  dealing  with 
these  findings,  we  need  to  look  in  more  detail  at  the  inception  scaling 
effects  that  had  already  been  reported.  Of  these,  certainly  the  work  by 
Parkin  and  Holl  (1953)  has  been  most  influential  in  providing  a  source  of 
information  used  by  co-workers  in  formulating  theories.  There  is  a 
pronounced  Reynolds-number  effect  on  these  bodies,  and  a  "size"  effect  (or 
unit  Reynolds-number  effect)  showing  that  additional  parameters  are 
important  (Figure  3).  An  even  larger  size  effect  was  produced  in  an 
inception  study  on  hydrofoils,  Parkin  (1952),  as  seen  in  Figure  4.  These 
effects  were  so  large  that  it  is  easy  to  understand  the  intense  effort  to 
explain  them  by  many  laboratories  and  investigators.  Despite  many  attempts, 
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Figure  3  -  Cavitation  Index  on  Hemispherical  and  1.5-Cal  Ogive  Bodies  as  a  Function  of  Reynolds 
Number  (From  Parkin  and  Holl  1953:  The  point  of  disappearance  of  the  cavitation  is  used 
as  the  index  in  this  case  and  this  is  termed  "desinent"  cavitation  by  Holl.) 


Cavitation  Index  on  a  Series  of  Geometrically  Similar  Hydrofoils  Showing  a  Size  Effect 

(Parkin  1952) 


no  simple  correlations  with  surface  tension  or  "time  of  travel"  have 
proved  successful.  And,  as  suggested,  the  subsequent  theories  based 
largely  on  Parkin  and  Kermeen's  1953  observations  have  had  difficulties. 

As  the  final  example  of  Table  1,  we  may  cite  again,  that  of  Holl  and 
Kornhauser  (1970).  They,  like  van  der  Walle  (1962),  used  the  conceptual 
idea  of  nuclei  stabilization  on  the  body  by  pressure  gradient  in  a  region 
of  tension,  thereby,  permitting  growth  through  the  Plesset-Zwick  bubble 
mechanics  to  explain  the  thermal  effects  experienced  in  hot  water  and  some 
non-aqueous  fluids.  Unlike  many  coworkers,  they  carried  out  many  experi¬ 
ments,  only  to  find  a  negative  correlation  for  the  proposed  mechanism. 

All  of  these  theories,  and  others,  are  summarized  by  Holl  in  his 
exhaustive  review  article  in  "Cavitation  State  of  Knowledge,"  (Robertson 
and  Wislicenus  1969).  At  that  time,  the  International  Towing  Tank 
Conference  (ITTC)  round  robin  tests  were,  regrettably,  not  very  widely 
known. 

The  ITTC  Tests 

Under  the  auspices  of  the  ITTC,  tests  on  a  flat-faced  ellipsoidal  body 
were  carried  out  in  many  different  laboratories  throughout  the  world,  with 
the  apparently  chaotic  results  seen  in  Figure  5.  There  was  no  doubt  that 
these  tests  were  extremely  thought-provoking,  even  if  they  did  show  a 
potentially  embarrassing  lack  of  uniformity  in  a  standard  test.  Some  of 
the  oddities  of  these  results,  namely  »  -c  were  shown  by  Holl  to 

be  due  to  attached  cavities  in  more-or-less  equilibrium  with  the  dissolved 

air  in  the  tunnel  water.  Then  p  .  =  p  where  p  _  is  the  gas  satu- 

min  sat  *sat  & 

ration  pressure,  so  that  a.  -  p  /q  and  as  q  decreases  the  cavitation 

l  sat 

index  may  become  quite  large.  At  about  this  time  it  was  suggested  by 
Johnsson  (1969)  that  perhaps  these  differences  were  due  to  a  laminar 
separation.  Peterson  (1969),  had,  in  fact,  already  observed  a  laminar 
separation  on  a  similar  type  of  ellipsoidal  body,  but  at  tunnel  speeds 
less  than  those  of  his  subsequent  inception  tests.  A  most  graphic  and 
telling  comparison  was  made  by  Johnsson  (ibid)  who  showed  a  photograph  of 
the  cavitation  actually  seen  on  some  of  the  ITTC  test  bodies.  Then  part 
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of  the  problem  became  immediately  clear.  We  see  in  Figure  6  a  reproduction 
of  this  photograph  in  which  it  can  be  seen  that  quite  different  kinds  of 
cavitation  are  present  on  those  test  bodies  and  a  closeup  photograph  taken 
in  the  High  Speed  Water  Tunnel  (HSWT)  at  Caltech  (Figure  7),  reveals  an 
attached,  smooth  form  of  cavitation  with  no  travelling  bubbles  whatever! 

Forms  of  Cavitation 

The  photograph  of  Figure  7  is  certainly  far  removed  from  the  concept 
of  travelling-bubble  cavitation,  and  although  it  is  quite  different  from 
many  contemporary  photographs  in  other  facilities,  it  is  by  no  means 
unique.  Eventually,  by  observation  on  this  and  other  bodies,  such  as  the 
hydrofoil  of  Figure  8,  we,  together  with  other  workers,  became  convinced 
that  three  more-or-less  different  kinds  of  cavitation  could  usually  be 
distinguished.  These  are:* 

1.  Travelling-bubble  cavitation. 

2.  Sheet  or  band  cavitation. 

3.  Attached  spots,  "fingers”  or  "wedges"  of 
cavitation. 

All  of  these  forms  can  be  seen  in  Figure  6.  In  the  light  of  all  these 
observations  and  suggestions  it  seemed  an  appropriate  time  to  consider,  in 
more  detail,  the  actual  viscous  flow  about  these  various  bodies  that  had 
previously  been  studied. 

VISCOUS  EFFECTS  ON  CAVITATION  INCEPTION 
Experimental  Methods 

As  in  aerodynamic  flows,  we  expect,  in  hydrodynamic  applications, 
changes  to  be  brought  about  by  viscosity  in  surface  boundary  layers.  We 
expect  to  be  wary  of  laminar  separation  in  adverse  pressure  gradients  and 
to  anticipate  laminar-to-turbulent-f low  transition  of  attached  boundary 
layers  at  sufficiently  high  Reynolds  numbers.  Also  we  know,  in  a  general 
way,  that  although  transition  is  a  complex  process,  a  precise  knowledge  of 
its  location  is  rarely  required  except  for  drag  calculations.  But  to  know, 

*Vortex  cavitation  also  occurs  on  lifting  surfaces  but  is  not  included 
in  the  present  discussion. 


in  more  detail,  which  of  these  features  (as  well  as  others  not  mentioned)  is 
actually  present  in  particular  hydrodynamic  experiments  requires  some  means 
of  boundary-layer  flow  visualization.  Surface  films  of  various  types  have 
traditionally  been  used  in  aerodynamic  applications  for  this  purpose  and 
in  hydrodynamic  work  oil  films  have  been  found  useful  for  detection  of 
laminar  separation  and  transition  (see,  e.g.,  the  review  Acosta  and  Parkin 
1975).  But  a  much  more  sensitive  method  for  water  is  schlieren  photogra¬ 
phy.  It  is,  of  course,  necessary  in  this  method  to  provide  a  density 
contrast;  this  can  be  achieved  by  injection  of  another  liquid  or,  more 
simply,  by  heating  or  cooling  the  surface  layers.  One  sees  then,  in  the 
usual  schlieren  setup  the  thermal  boundary  layer  which  is  for  water  at 
room  temperature  about  one-half  the  mechanical  boundary  layer'  thickness. 

We  were  concerned  at  first  about  the  effects  of  this  heating  or  cool¬ 
ing  on  the  boundary  layers  themselves.  Fortunately,  this  problem  had 
already  been  exhaustively  treated  by  Wazzan  et  al.  (1968),  (see  also  Wazzan 
and  Gazley  1977)  to  show  that  under  the  expected  laboratory  conditions  the 
effect  was  negligible  near  a  laminar  separation  but  that  some  delay  in  a 
turbulent  transition  might  be  expected.  (We  shall  return  to  this  point 
later.)  A  schematic  arrangement  is  shown  in  Figure  9;  this  is  basically 
the  same  as  that  of  Arakeri  et  al.  (1973)  but  with  an  improved  spark 
source.  A  photograph  of  the  experimental  arrangement  in  the  Caltech  Low 
Turbulence  Water  Tunnel  (LTWT)  may  be  seen  in  Figure  10. 

Laminar  Separation 

The  schlieren  technique  was  surprisingly  easy  to  apply.  For  example. 
Figure  11  shows  a  photograph  of  the  flow  past  a  2-in.  diameter  hemisphere 
body  at  a  flow  speed  of  about  22  ft/sec,  well  within  the  range  of  experi¬ 
mental  Reynolds  numbers  of  the  data  of  Figure  3.  The  clear  laminar 
separation  on  this  body  is  easily  seen,  together  with  the  turbulent  re¬ 
attachment  of  the  free  shear  layer  further  downstream.  The  results  really 
exceeded  our  expectations  for  it  could  be  quickly  confirmed  that  a  laminar 
separation  was  present  on  the  hemispherical  body  up  to  the  top  speed  of 
about  60  ft/sec  in  the  HSWT  at  Caltech  (a  body  Reynolds  number  of  about 
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Figure  9  -  Schlieren  Flow  Visualization  System 
(Lindgren  and  Johnsson  1966) 


(Lindgren  and  Johnsson  1966) 


10^)  and  that  laminar  separations  were  present  on  the  ITTC  test  body  and 
the  modified  ellipsoidal  head  form  previously  tested  by  Peterson  at  the 
David  W.  Taylor  Naval  Ship  Research  and  Development  Center  (see  Arakeri  and 
Acosta  1973,  1976). 

At  first,  it  was  surprising  to  realize  that  laminar  separation  could 
exist  on  these  well-known  test  bodies  to  comparatively  high  Reynolds 
numbers  because  of  contrary  assumptions  found  in  some  of  the  cavitation 
literature.  Yet,  even  the  simplest  approximate  laminar  boundary- layer 
calculations*  showed  that,  on  these  bodies,  a  separation  was  to  be  expect¬ 
ed.  It  then  seemed  evident  that  the  presence  or  absence  of  this  laminar 
separation  would  be  extremely  important  in  cavitation  inception. 

Inception  Observations  within  the  Laminar 
Separation 

With  the  schlieren  system  in  operation,  it  became  possible  to  carry 
out  simultaneous  observations  of  incention  and  the  thermal  boundary  layer. 
The  hemisphere  body  proved  to  be  a  good  candidate  for  this  purpose.  The 
development  of  cavitation  on  this  body  may  be  seen  in  Figure  12.  Cavi¬ 
tation  is  seen  to  commence  within  the  reattachment  region  of  the  laminar 
separation.  Many  small  bubbles  are  seen  within  the  recirculation  region 
of  the  separation;  these  are  presumably  the  same  small  bubbles  seen  in 
Parkin  and  Kermeen* s  photographs.  Motion  pictures  were  also  taken  of  the 
flow  in  this  region,  at  the  same  magnification  of  the  schlieren  photographs. 
From  these,  it  was  inescapable  that  these  small  bubbles  were  in  the  re¬ 
circulation  region  of  the  separation,  but  it  was  not  possible  to  establish  a 
direct  relation  between  these  very  small  bubbles  and  the  larger  macroscopic 
cavitation  patches  seen  further  downstream. 

With  continued  reduction  in  pressure,  the  macroscopic  region  of  cavi¬ 
tation  grows  to  eventually  result  in  the  smooth,  clear  viscous-cavitation 
separation  seen  in  the  last  of  Figure  12.  This  is  a  new  type  of  separation 
phenomenon,  somewhat  related  to  separation  in  lubrication.  (Arakeri 
(1975a)  was  able  to  make  semi-empirical  correlations  for  the  location  of 
this  separation.)  A  similar  sequence  of  events  was  found  for  cavitation 


*For  example  the  axisymmetric  version  of  Thwaites  method,  see  "Laminar 
Boundary  Layers,"  Oxford,  Rosenhead  (ed.)  pp.  430-432. 
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(b)  a  =  0.70 


Figure  12  -  The  Form  and  Extent  of  Cavitation  Originating  Within  the 
Viscous  Separated  Region  of  the  Hemispherical  Nose  at  Three 
Levels  of  Tunnel  Pressure  (The  freestream  speed  is 
AO  ft/ sec  and  the  body  Reynolds  number  is 

6  x  10^.  The  dark  patches  above  the 
model  outline  are  the  cavitating 
areas.  Arakeri  and  Acosta 
1973) 


development  on  the  ITTC  body.  This  body,  however,  has  a  much  smaller 
laminar-separation  region  and  in  the  Caltech  HSWT  the  progress  from 
initial  inception  to  the  macroscopic  cavitation  seen  in  Figure  7  is 
practically  instantaneous.  But,  again,  careful  study  shows  that  the  ter¬ 
minus  of  the  laminar  region  is  the  location  of  the  most  intense  cavitation, 
and  it  is  this  location  where  cavitation  is  seen  to  disappear  as  pressure 
is  raised. 

We  can  see  that  a  laminar-separation  region  has  an  important  role  in 
cavitation  inception.  At  the  very  least,  the  static  pressure  there  pro¬ 
vides  the  reference  pressure  at  which  inception  occurs.  It  would  then  seem 
most  plausible  that 


a 


i 


(6) 


or 


°i  =  "cp  (7) 

preattach 

i.e.,  we  identify  the  inception  index  with  the  pressure  coefficient  at 
separation  or  reattachment,  the  latter  being  Reynolds  number  dependent. 

In  fact,  this  suggestion  had  previously  been  made  by  Bailey  (1970)  based  on 
hydrofoil  tests.  Here  we  exhibit  tests  on  the  two  bodies,  hemisphere  and 
ITTC,  made  in  the  Caltech  HSWT  in  Figures  13  and  14,  respectively,  in  which 
it  may  be  seen  that  the  rule  is  better  than 


a . 

i 


-c 

p  . 
min 


(8) 


and  that  the  inception  index  approaches  -c  as  speed  increases.  But, 

Ps 

there  are  differences  in  the  two  sets  of  results:  namely, 
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TUNNEL  VELOCITY  -  U  IN  FEET  PER  SECOND 


(9) 


21  -cp  hemisphere  body 


and 


ai  £  -c  ITTC  body  (10) 

ps 

Although  there  is  some  improvement  when  the  reattachment  pressure  co¬ 
efficient  is  used,*  the  two  bodies  are  seen  to  behave  quite  differently. 

And  the  original  observations  of  Parkin  and  Kermeen  (namely  that  the  local 
pressure  exceeds  vapor  pressure  yet  the  cavitation  still  occurs)  on  the 
hemisphere  body  are  confirmed  once  again. 

Before  moving  on  to  questions  raised  by  these  data,  it  is  appropriate 
to  point  out  that  in  this  section  we  have  been  dealing  with  attached  forms 
of  cavitation  which  are  seen  to  occur  within  a  boundary  layer  separation. 

It  seems  clear  that  the  sheets,  or  band  forms  of  cavitation,  are  due  to  this 
separation.  In  the  Caltech  experiments  so  far  reported,  this  type  of  cavi¬ 
tation  occurs  without  travelling  bubble  cavitation.  But  we  see,  for  example 
on  the  ITTC  body  in  Figure  6,  that  other  forms  of  cavitation  also  occur  and 

it  should  not  be  supposed  that  the  approximate  correlation,  =  -c  ,  is 

Ps 

valid  for  them.  Indeed,  the  great  variability  of  the  ITTC  inception 
results  of  Figure  5  makes  it  plain  that  each  of  these  forms  of  cavitation 
inception  will  have  its  own  governing  dynamics.  Of  course,  one  is  left 
with  the  problem  of  determining,  in  a  particular  flow  situation,  the 
particular  form  of  cavitation,  so  that  the  appropriate  physical  model  may 
be  used  to  determine  the  scaling  procedure. 

Conditions  within  the  Separation 

There  is  still  the  problem  of  explaining  cavitation  on  the  hemisphere 
body  when  the  local  pressure  (in  the  separation  region)  exceeds  vapor 
pressure.  (We  exclude  gaseous  cavitation  here  by  the  slowness  of  the 
diffusion  process.)  The  implication  is  clear  that  there  must  be  local. 


*See  Arakeri  (1973)  for  these  measurements. 
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transient  pressure  fluctuations  sufficiently  strong  to  create  liquid 
tensions  sufficient  for  microbubble  growth.  From  the  data  of  Figure  13, 
for  example,  a  fluctuation  magnitude  of  about  ten  percent  of  the  dynamic 
pressure  would  be  required  to  establish  cavitation  on  the  hemisphere  body. 
These  fluctuations  would  be  expected  to  occur  at  the  reattachment  of  the 
separation  region,  as  can  readily  be  appreciated  from  Figure  15.  Rather 
surprisingly,  the  literature  is  very  scarce  on  this  feature  of  real  fluid 
flows.  It  was  possible,  however,  to  install  a  single  point  pressure 
transducer  on  the  hemisphere  body  of  Figure  15  (Arakeri  1975b).  The 
results  were  most  interesting  and  show  that  rms  and  peak  fluctuations  are 
certainly  large  enough  in  themselves  to  cause  liquid  tensions  in  the 
regions  of  observed  inception.  Later,  these  findings  were  confirmed,  with 
even  more  dramatic  and  precise  results  on  a  larger  body,  by  Huang  and 
Hannan  (1975).  Some  of  their  findings  are  shown  in  Figure  16,  to  demon¬ 
strate  that  the  fluctuating  pressure  environment  is  sufficiently  strong  to 
make  it  very  likely  for  cavitation  to  occur  on  these  bodies  at  local  pres¬ 
sures  slightly  exceeding  the  liquid  vapor  pressure. 

From  these  findings  then,  it  is  tempting  to  account  for  these  tran¬ 
sient  fluctuations  by  putting  the  inception  index  in  the  form 


where  c  is  a  local  pressure  coefficient  where  inception  is  observed, 
Ploc 

and  c  is  the  fluctuating  pressure-coefficient  amplitude  there.  This  type 
Pt 

of  representation  indeed  makes  the  hemisphere  data  much  more  plausible  be¬ 
cause  the  c  ,  measured  near  the  end  of  the  laminar  separation  bubble, 
ploc 

(and  perhaps  the  fluctuating  component  t oo )  is  strongly  Reynolds  number 

sensitive  (Huang  and  Hannan  ibid) .  This  prescription  ignores  the  physical 

state  of  the  fluid,  however,  so  that  some  caution  should  be  exercised  in 

seeking  correlations  based  solely  on  Equation  (11),  and  does  not  explain 

the  contrary  finding  of  Figure  14,  namely  that  a ^  <  -c  or  -c 

Ps  Preattach 
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Nevertheless,  there  is  considerable  appeal  of  the  idea  that  a  more 
realistic  appraisal  of  the  actual  features  of  the  real  viscous  flow  around 
bodies  of  interest  may  help  considerably  in  improving  our  knowledge  of 
cavitation  inception  and  the  development  of  better  scaling  laws.  This,  in 
turn,  leads  us  to  reflect  on  just  what  the  locations  of  inception  might  be 


(and  hence  the  value  of  c 


there)  and  the  strength  of  any  unsteady 


fluctuations  there — without,  for  the  moment,  worrying  about  how  the  cavi¬ 
tation  might  originate.  With  this  in  mind,  we  summarize  some  of  these  real 
fluid  features  that  might  influence  inception  on  smooth  bodies: 

1.  The  laminar-separation  reattachment  region. 

2.  A  "natural"  transition. 

3.  Stimulated  or  "tripped"  boundary  layers. 

4.  Freestream  effects,  e.g.,  turbulence,  particulates, 
polymer  solutions. 

We  cannot  deal  extensively  with  these  well-known  fluid  mechanics 
problems;  much  of  the  presently- known  structure  of  laminar-separation 
bubbles  may  be  found  in  Gaster  (1966)  and  the  principal  new  features  here 
are  the  astoundingly  large  pressure  fluctuations  as  already  mentioned. 

There  has  been  recent  progress  in  determining  the  length  of  these 
separations,  and  of  the  turbulent  transition  process  seen  on  the  free  shear 
layer  of  these  flows  (Van  Ingen  1975).  But  it  is  a  well-known  fact  that 
the  steady,  laminar  flow  on  smooth  bodies  gives  way  to  a  turbulent  flow  at 
a  sufficiently  large  Reynolds  number  and  that,  if  this  transition  is 
completed  upstream  of  a  laminar  separation,  this  separation  will  disappear. 
Flow  transition  may  be  stimulated  by  many  methods  if  the  Reynolds  number  is 
not  high  enough  for  a  "natural"  one — as  it,  indeed,  has  been  done  histori¬ 
cally  in  naval  hydrodynamics.  And,  if  by  this  process,  a  prior  laminar 
separation  is  forstalled,  one  would  expect  major  consequences  for  the  cavi¬ 
tation  inception  process  there.  Other  traditional,  and  some  nontraditional , 
issues  naturally  arise  when  transition  in  laboratory  facilities  occurs. 

Of  these,  surely  that  of  freestream  turbulence  has  been  one  of  the  most  re¬ 
curring  ones.  In  hydrodynamic  applications,  solute  molecules  and  particu¬ 
late  suspensions  may  have  profound  effects  on  both  laminar  and  turbulent 
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flow  structure.*  Of  these,  the  long  chain  'polyox'  polymer  has  received 
much  attention  in  the  hydrodynamic  community  for  its  astounding  effect  on 
turbulent  drag  and  heat  transfer  reduction.  It  would  be  surprising,  in  a 
way,  if  there  were  not  an  equivalent  effect  for  cavitation.  We  now  take  up 
some  of  the  features  briefly,  but  with  primary  reference  to  cavitation. 

Transition  and  Cavitation 

Here  we  need  to  have  some  idea  of  the  location  of  "natural"  transition 
of  attached  boundary  layers  to  determine  the  "local"  pressure  needed  for  a 
cavitation- inception  measurement.  This  is  a  basic  "reference"  problem  in 
fluid  mechanics  that  does  not  appear  yet  to  be  in  a  satisfactory  state 
(Mack  1978) .  At  best,  there  appear  to  be  semi- empirical  methods  based  upon 
the  use  of  linear  stability  theory.  In  its  present  form,  the  concept,  as 
explained  in  the  comprehensive  report  by  Wazzan,  Okamura,  and  Smith  (1968), 
is  to  calculate  the  factor  by  which  certain  disturbances  in  the  laminar 
boundary  layer  may  be  amplified.  As  an  empirical  observation,  transition 
may  be  said  to  occur  when  this  factor  reaches  a  certain  level.  The 
difficulties  and  objections  to  this  approach  are  reviewed  at  length  by 
Reshotko  (1976)  and  later  by  Mack  (ibid)  and  Van  Ingen  (ibid).  These 
center  on  the  fact  that  the  ambient  unsteady  disturbances  in  the  flow  are 
largely  unknown,  and  that  the  "receptiveness"  of  the  boundary  layer  to 
these  disturbances  is  not  known  either,  a  concept  due  to  Morkovin  (see 
Reshotko  1976).  It  does  not  necessarily  follow  that  transition  is  due  only 
to  growing,  unstable  two-dimensional  waves;  this  process  may  be  circumvent¬ 
ed  by  a  roughness  element,  for  example.  Nevertheless,  this  method,  the  e11 
method,  as  it  is  called,  has  value  in  comparing  different  flows  in  a  fixed 
environment,  and  is  in  any  case  the  most  quantitative  method  now  available. 
For  example,  simplified  prediction  schemes  based  on  these  stability  methods 
can  now  readily  account  for  transition  of  heated  water  boundary  layers 
(Wazzan  and  Gazley  1978) . 

With  the  linear  stability  method  of  Smith,  Wazzan,  and  coworkers,  it 
is  readily  possible  to  calculate  the  transition  location  (or  the  start  of 

*Suspended  'dust'  particles  may  have  some  similar  effects  in  gas  flows. 
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transition)  for  a  range  of  amplification  levels.  Values  of  the  level 

7  13  9 

chosen  range  from  e  (Arakeri  1975b)  to  e  (Huang  and  Hannan  1975)  with  e 

as  the  original  value  chosen  by  Smith.  Many  such  transition  estimates  have 
now  been  made,  e.g.,  on  bodies  of  interest  in  naval  hydrodynamics  (Kaups 
1974  and  Power  1977),  but  Arakeri  (1973)  was  evidently  among  the  first  to 
apply  this  idea  to  cavitation  inception.  Calculations  were  made  of  the 
pressure  coefficient  at  transition  (using  the  amplification  level  of  e^), 
and  the  results  so  obtained  were  compared  with  the  desinent  cavitation  ob¬ 
servations  of  Parkin  and  Holl  (1953)  on  a  1.5-cal  ogive.  These  results, 
reproduced  in  Figure  17,  are  extremely  suggestive.  Later,  it  was  possible 
to  measure  surface  pressure  fluctuations  at  a  point  on  this  same  body 
(Figure  18)  to  reveal  large  peak  pressure  fluctuations  at  the  site  of 
transition.  The  spectra  of  these  same  fluctuations  were  compared  with 
the  expectations  of  the  linear  stability  calculations  (Figure  19)  to 
show,  remarkably,  that  the  most  sensitive  frequencies  according  to  the 
theory  were  actually  observed.  So  it  was  with  some  anticipation  that 
spark  schlieren  photographs  were  taken  of  this  body  (Figure  20) .  These 
appeared  to  reveal  cavitation  inception  at  the  observed  position  of  what 
are  interpreted  to  be  boundary  layer  transition  disturbances  (Arakeri  and 
Acosta  1974). 

More  recent  work  of  McCarthy  et  al.  (1976)  has  shown  a  similar 
favorable  agreement  of  surface  (hot  film)  fluctuation  frequency  response 
with  the  estimates  made  from  linear  stability  calculations.  Perhaps  then 
the  plausibility  of  cavitation  inception  correlations  with  the  location  of 
transition,  i.e.,  a  relation  such  as  Equation  (11),  is  heightened  with  this 
closer  association  of  transition  as  a  site  of  fairly  large  disturbances. 
Direct  confirmation  of  transition-connected  cavitation  on  small  bodies 
having  deeper  minimum  pressures  present  some  experimental  problems  because 
of  the  small  range  of  speed  and  pressure  variables  available  for  cavitation 
index  and  Reynolds  number  excursions.  One  such  body  that  has  proved  useful 
in  laboratory  comparisons  is  the  modified  ellipsoidal  flat  faced  head  form 
similar  to  the  ITTC  body  except  with  a  deeper  pressure  minimum.  It  was 
used  in  cavitation  studies  first  by  Peterson  (1969)  and  then  later  by 
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Pressure  Fluctuations  on  the  Surface  of  a  1.5~Cal  Ogive  at  the  Arc-Diameter  Ratio  of 

(Arakeri  1975b) 


Figure  19  -  Comparison  of  Measured  Frequency  Distribution  with  that  Calculated  from  Linear  Stability 

Results  (Arakeri  1975b) 


Brockett  (1972).  This  body,  termed  the  NSRDC  body,  like  the  ITTC  body, 
possesses  a  laminar  separation  up  to  a  critical  Reynolds  number  based  on  a 
diameter  of  about  0.46  x  10  at  which  both  the  e°  method  (n  =  7)  and 
observations  show  transition  to  occur,  thereby  forestalling  the  presence 
of  separation. 

For  Reynolds  numbers  less  than  this  critical  one  the  appearance  of 
cavitation  on  this  body  is  very  similar  to  that  on  the  ITTC  body,  i.e.,  a 
smooth  band  cavity  is  seen  (in  the  Caltech  HSWT) .  At  greater  Reynolds 
numbers  this  band  type  gives  way  to  patchy  regions  of  foamy  cavitation 
bubbles,  intermittently  spaced  around  the  periphery  with  the  smooth  sheet 

cavity.  Finally,  at  greater  Reynolds  numbers  of  say  7  x  10^,  spots  or 

wedges  of  cavitation  are  seen.  The  inception  index  on  one  such  body  has 
been  determined  in  two  different  facilities  with  quite  different  results 
(Arakeri  and  Acosta  1976).  In  the  Caltech  facility  the  correlation 

-  c  was  found  to  be  valid  for  Reynolds  numbers  less  than  the  critical 

one.  A  natural  transition,  precluding  the  separation,  existed  for  Reynolds 
numbers  greater  than  about  5  x  10^.  The  observed  inception  index  remained 
quite  close  to  the  pressure  coefficient  at  the  site  of  the  predicted 
transition  location  (Figure  21)  though  there  was  hardly  any  change  in  the 
pressure  coefficient  from  the  separation  value. 

Spot  or  wedge  cavitation  was  also  observed  to  occur  at  the  greatest 

values  of  Reynolds  number  originating  from  positions  near  c  .  These 

Pmin 

indices  are  quite  different  from  the  attached  cavitation  associated  with 
the  presumed  location  of  transition. 

However,  more  recent  studies,  by  Huang  and  Santelli  (1977),  show  that 
on  one  body  not  subject  to  a  laminar  separation,  travelling-bubble  inception 
was,  in  fact,  observed  at  the  estimated  position  of  transition.  From  these 
observations,  they  deduce  an  inception  index  behavior  of  the  form  of 
Equation  (11)  with  the  local  pressure  coefficient  being  that  of  transition 
as  in  Figure  17  and  with  the  unsteady  term  being  about  -0.06.  This  latter 
is,  however,  air  content  sensitive,  which  is  a  subject  we  shall  return  to 
later. 
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-  Inception  and  Desinent  Cavitation  Index  for  the  NSRDC  Body  in  the  Caltech  HSWT  at  an 
r  Content  of  10.4  ppm  (The  numbers  in  parentheses  represent  the  number  of  data 
>oints  and  the  height  of  the  bars,  the  range  of  the  data.  Laminar  separation 
was  not  observed  greater  than  about  28  ft/sec.)  (Arakeri  and  Acosta  1976) 


Stimulated  Boundary  Layers 

From  the  previous  discussion,  we  would  expect  that  any  effect  on  cavi¬ 
tation  of  a  laminar  separation  would  be  eliminated  if  the  laminar  boundary 
is  sufficiently  "stimulated"  to  transition.  This  stimulation  may  be  via  a 
variety  of  freestream  factors,  but  in  the  laboratory  it  is  easiest  to  use 
a  mechanical  "isolated  roughness"  or  boundary  layer  "trip."  Studies  were 
carried  out,  first  on  the  hemisphere  body  (Arakeri  and  Acosta  1976),  with 
trips  of  tape  and  then  with  machined  steps.  These  trips  were  rather  large 
(several  times  larger  than  the  boundary  layer  displacement  thickness)  and 
were  located  near  the  nose  to  forestall  cavitation  at  the  trip.  At  low 
speeds  a  laminar  separation  was  still  observed  (Figure  22),  but  at  higher 
ones,  about  30  ft/sec,  the  separation  gave  way  as  transition  occurred.  The 
results  of  this  tripping  were  indeed  spectacular,  as  the  normal  course  of  a 
cavitation  experiment  was  completely  reversed.  We  see  in  Figure  23  a 
sequence  of  photographs  of  cavitation  on  a  hemisphere  body  taken  as  speed 
is  increased.  As  this  happens  the  established  laminar  separation  disappear 
and  with  it,  in  the  Caltech  HSWT,  the  cavitation,  tool  In  fact,  cavitation 
inception  indices  show  a  decrease  with  tunnel  speed  (Figure  24)  instead  of 
the  usual  trend.  We  should  emphasize  that  the  form  of  cavitation  in  this 
tunnel  was  always  that  of  an  attached  cavity;  freestream  travelling  cavi¬ 
tation  events  were  extremely  rare.  However,  this  kind  of  cavitation 
observation  is  not  unique  to  the  Caltech  facility,  as  the  interesting 
photograph  of  Figure  25  shows.  More  recently,  boundary-layer  stimulation 
has  been  proposed  as  a  standard  laboratory  technique  in  cavitation  testing 
to  avoid  unwanted  scale  effects  from  laminar  separation  (Kuiper  1978).  As 
the  findings  of  Figure  23  shows  there  still  may  be  other  important  factors. 

These  findings  on  tripped  boundary  layers  raise  many  issues  concerning 
inception  scaling  laws.  These  issues  are  of  two  types;  in  the  first,  we 
have  (as  already  sketched)  all  of  those  questions  affecting  the  fully- 
wetted  viscous  flow  itself.  The  second  type  may  be  said  to  be  those  having 
a  direct  bearing  on  the  cavitation  process  itself.  Among  these  is  the 
presence  (or  absence)  of  nuclei,  and  the  state  of  the  liquid  (i.e.,  thermo¬ 
dynamic  factors).  There  is  also  the  large  effect  of  polymer  solutions  on 


Figure  23  -  A  Two-Inch  Hemisphere  Body  with  Trips  Showing  the 
Disappearance  of  Cavitation  as  Tunnel  Speed  is 
Increased  and  the  Cavitation  Index  Decreased 
(The  left  side  sequence  has  a  trip  height 
of  0.012  inches  and  the  right  side 
has  a  trip  of  0.005  inches, 

Arakeri  and  Acosta  1976) 
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portion  of  the  body  the  normal  band  cavitation  can  be 
Photo  courtesy  of  Professor  J.  W.  Holl,  Penn  State 
University. ) 


decreasing,  perhaps  by  one-half,  the  cavitation  inception  index  in  jets  and 
the  flow  over  bodies  reported  by  many  workers  (e.g.,  Ellis  and  coworkers, 
Holl,  Arndt  and  coworkers,  Hoyt,  van  der  Meulen;  see  van  der  Meulen  1974, 
Gates  1977  for  review).  Despite  the  initial  conclusions  of  Ellis  (Ellis 
et  al.  1970)  that  bubble  mechanics  were  not  a  factor  in  this  suppression, 
considerable  uncertainty  as  to  the  basic  cause  of  the  effect  remained. 

The  Polymer  Effect  in  Cavitation 

It  is  certainly  fair  to  say  that  the  presence  of  long-chain  polymers 
in  dilute  liquid  solutions  has  had  a  vitalizing  effect  on  liquid  fluid 
mechanics.  The  many  reasons  for  this  development  are  well  covered  by  the 
extensive  review  of  Hoyt  (1972).  But,  as  sketched  in  the  foregoing  dis¬ 
cussion,  there  is  an  additional  unique  effect  of  polymer  solutions  on  cavi¬ 
tation  inception.  Two  kinds  of  flows  have  been  extensively  pursued, 
cavitation  inception  on  bodies  (as  discussed  herein)  and  in  jets.  The 
magnitude  of  the  effect  appears  to  be  equivalent;  we  discuss  only  smooth 
body  inception  in  the  following  sections. 

It  seems  most  natural  in  discussing  these  types  of  fluid  to  note  the 
appearance  of  particular  phenomena  with  a  characteristic  time  ratio,  such 
as  the  ratio  of  the  time  required  for  the  fluid  to  move  a  certain  length, 
to  a  characteristic  polymer  molecular  relaxation  time.  In  fact  such  a 
relation  was  used  by  Arndt  (Arndt  et  al.  1976)  to  correlate  such  suppression 
results,  and  it  appeared  to  work  best  if  a  relevant  boundary  layer  length 
was  used  as  the  length  parameter.  At  about  the  time  of  the  17th  ATTC 
conference  (Morgan  1974)  it  was  suggested  by  several  workers  including  J.W. 
Holl  that,  based  on  the  appearance  of  the  cavitation  band  at  inception, 
the  polymer  seemed  to  cause  the  interface  to  become  turbulent  with  the 
implication  that  transition  had  occurred.  But  this  demonstration  remained 
for  van  der  Meulen  (1976)  to  carry  out.  He  showed,  by  careful  and  ingeni¬ 
ous  adaptation  of  a  schlieren-holographic  technique,  that  the  polymer  fluid 
appears  to  trigger  a  turbulent  transition  and  that,  on  the  hemispherical 
body  (for  example),  this  transition  removed  the  laminar  separation.  Then 
one  would  have,  presumably,  a  situation  rather  similar  to  that  of  a  body 
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normally  having  a  laminar  separation,  but  with  a  tripped  boundary  layer, 
and  equivalent  effects  might  therefore  be  anticipated.  Thus,  as  in 
Figures  23  and  24,  a  cavitation  suppression  effect  should  be  expected  from 
the  polymer  fluid  if  the  main  effect  of  the  polymer  is  on  the  boundary 
layer.  Indeed,  the  magnitude  of  van  der  Meulen's  suppression  effect  are 
similar  to,  but  not  as  great  as,  those  of  Figure  24. 

van  der  Meulen  further  reasoned  that  a  body  not  normally  possessing 
a  separation  should  not  exhibit  a  large  suppression  effect.  To  this  end, 
he  adopted  a  body  shape  previously  suggested  for  hydrodynamic  cavitation 
research  by  Schiebe  (Schiebe  1972);  namely,  a  body  having  nearly  the  same 


minimum  pressure  coefficient  as  the  hemisphere  (c 


=  -0.75)  but  without 


^min 

a  predicted  laminar  separation,  van  der  Meulen's  experiments  showed  little 
effect  of  the  polymer  injectant  on  the  cavitation- inception  behavior  of 
this  body.  Furthermore,  his  holographic  schlieren  observations  confirmed 
the  predicted  absence  of  separation.  He  was  also  able  to  show  that  the 
position  of  transition  estimated  from  these  observations  agreed  on  the 
whole  quite  well  with  the  linear  stability  calculations.  Interestingly 
in  his  test  facility,  travelling-bubble  cavitation  was  the  predominant 
form  of  cavitation  observed  on  this  Schiebe  body.  In  a  followup  report 
(van  der  Meulen  1978)  he  summarizes  his  opinions  about  the  hemisphere  and 
Schiebe  bodies  thusly,  "Since  the  influence  of  polymer  additives  is  to 
suppress  laminar  boundary- layer  separation  on  the  hemispherical  nose,  the 
strong  pressure  fluctuations  occurring  at  the  position  of  transition  and 
reattachment  of  the  separated  shear  layers;  and  being  the  principal 
mechanism  for  cavitation  inception,  are  eliminated  and  cavitation  inception 
will  start  at  much  lower  pressures." 

Much  the  same  conclusions  had  been  reached  at  about  the  same  time  by 
Gates  (Gates  1977,  Gates  et  al.  1978)  on  the  mechanism  of  polymer  effect. 
Both  Gates  and  van  der  Meulen  injected  polymer  solutions  at  the  stagnation 
point.  Gates  was  able  to  observe  a  significant  polymer  effect  at  extremely 
low  injectant  rates  on  both  the  hemispherical  body  (Figure  26)  and  the 
Schiebe  body  (Figure  27).  In  these  photographs,  the  injectant  rate  is 
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Figure  26  -  Flow  Past  the  Hemisphere  Body  with  Injection  of  100  wppm 
Polyox  at  a  Reynolds  Number  of  3.9  x  10^  (The  dimensionless 
injectant  rates  are,  respectively:  (a)  no  injection, 

(b)  0.5  x  lo"6,  (c)  1.1  x  io-6,  (d)  1.7  x  lo~6,  and 
(e)  2.9  x  10  ,  Gates  and  Acosta  1978) 


Figure  27  -  Flow  Past  the  Schiebe  Body  at  a  Reynolds  Number  of  4.2  x  10 
with  Injection  of  500  wppm  Polyox  (The  dimensionless  injection 
parameters  are,  respectively:  (a)  zero,  (b)  2.3  x  10~6, 

(c)  1.5  x  10  ■*,  and  (d)  2.9  x  10  The  photographs 
are  each  0.2  body  diameters  in  length  centered  on 
arc  length  ratios  of  0.82,  0.75,  0.6,  0.53, 
respectively, Gates  and  Acosta  1978) 


(a) 

1 

■  vi . 
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normalized  to  give  the  weight  fraction  within  the  boundary  layer  displace¬ 
ment  thickness.  It  is  evident  that  the  very  small  concentration  of  only 
2.9  ppm  is  sufficient  to  remove  the  laminar  separation  on  the  hemisphere, 
and  significant  forward  movement  of  transition  on  the  Schiebe  body  is 
accomplished,  with  somewhat  larger  concentrations.  The  main  thrust  of  van 
der  Meulen's  findings,  then,  appear  to  be  confirmed,  and  the  cavitation  in¬ 
ception  inhibition  on  the  hemisphere  is  traced  to  an  overall  viscous  effect 

Turbulence  Levels 

We  now,  of  necessity,  touch  briefly  upon  this  ubiquitous  standard 
problem  in  experimental  fluid  mechanics;  namely,  the  influence  of  free- 
stream  turbulence  on  flow  past  bodies.  But,  with  the  restricted  emphasis 
herein  on  cavitation  inception,  our  concern  is  primarily  to  determine  if 
the  stimulated  transition  causes  a  laminar  separation  to  disappear  or  if 
there  is  any  direct  effect  on  cavitation  inception  of  bodies  not  having 
such  a  separation.  There  are,  of  course,  semi-empirical  rules  for  the 
turbulence  effect  (e.g..  Hall  and  Gibbings  1972,  Van  Driest  and  Blumer 
1963)  and  the  turbulence-modified  amplification  method  of  Mack  (1978).  But 
as  Reshotko  (1976)  has  made  clear  in  his  review,  the  coupling  between  these 
freestream  disturbances  and  boundary  layer  developments  is  unclear — 
particularly  for  large  disturbances  that  may  bypass  the  linear  growth  pro¬ 
cess.  Thus,  as  in  many  areas  of  naval  hydrodynamics,  recourse  must  still 
be  made  to  the  laboratory  to  study  particular  experimental  situations. 

It  has  not  been  customary  in  much  hydrodynamic  work  to  quantify 
turbulence  levels — let  alone  their  power  spectra — although  a  few  such  data 
on  several  water  tunnels  are  now  becoming  available  through  the  ITTC 
reports  (12th  ITTC  Cavitation  Committee).  Thus,  it  is  no  surprise  that 
there  is  no  account  of  the  effect  of  freestream  turbulence  on  flow  past 
bodies  of  interest  in  naval  hydrodynamic  applications.  One  problem  in 
large  hydrodynamic  facilities  has  been  the  expense  of  the  control  of  this 
real  fluid  feature.  For  this  purpose  a  water  tunnel  having  a  low  level  of 
freestream  turbulence  to  begin  with  is  essential.  The  key  features  of  one 
such  research  facility  are  shown  in  Figure  28.  Some  features,  particularly 
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convenient  for  this  kind  of  work,  are  the  large  (16:1)  contraction  ratio, 
the  vaned  elbows  on  each  leg,  the  ready  insertion  of  honeycomb  straight- 
eners  and  damping  screens  and,  finally,  the  insertion  of  turbulence¬ 
generating  grids  at  the  entrance  of  the  working  section  (see  Gates  1973 
for  a  full  description).  It  was  possible,  thereby,  to  vary  the  freestream 
level  by  nearly  a  factor  of  100  from  a  low  value  of  about  0.05  percent. 
(Unfortunately,  it  was  not  possible  to  carry  out  cavitation  experiments 
over  the  same  range.) 

It  was  anticipated  that  these  turbulence  levels  would  have  a  profound 
effect  on  bodies  near  transition  but  normally  having  a  separation,  such  as 
the  DTNSRDC  body.  We  see,  in  Figure  29,  that  a  turbulence  level  of  about 
one  percent  is  sufficient  for  transition  to  occur  upstream  of  the 
separation  on  this  body.  A  similar  upstream  progression  of  the  site  of 
transition  at  fixed  body  Reynolds  number  with  increasing  turbulence  level 
was  also  observed  on  the  Schiebe  body  (not  shown).  Indeed,  the  observa¬ 
tions  of  transition  location  together  with  those  of  van  der  Meulen  agree 
fairly  well  with  the  amplification  computations  of  Wazzan  and  Gazley,  Fig¬ 
ure  30.  Similar  observations  on  the  hemisphere  body  were  surprising  how¬ 
ever,  as  the  laminar  separation  appeared  to  be  totally  immune  to  the  turbu¬ 
lence  level!  This  unexpected  result  is  traced  (Gates  1977)  to  significant 
differences  in  the  sensitive  critical  frequency  of  the  respective  boundary 
layers  on  these  bodies  (deduced  from  linear  stability  theory),  the  highest 
frequency  being  for  the  hemisphere  body,  DTNSRDC  next,  and  then  Schiebe 
body.  Because  of  this  difference  in  frequency  and  reasoning  from  typical 
grid  turbulence  spectra,  it  is  argued  that  there  is  nearly  two  orders  of 
magnitude  more  disturbance  energy  available  for  boundary  stimulation  on 
the  DTNSRDC  body  than  on  the  hemisphere.  This  seems  a  plausible  explana¬ 
tion  and  shows  the  desirability  of  providing  freestream  spectral  data  in 
addition  to  intensity  levels  in  future  hydrodynamic  work. 
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Figure  29  -  The  Effect  of  Freestream  Turbulence  on  the  Flow  Past  the  NSRDC  Body 

(The  flow  is  from  right  to  left  at  a  body  Reynolds  number  of  1.6  *  10^.  The 
turbulence  intensity  is:  (a)  0.05  percent,  (b)  0.65,  (c)  1.1,  (d)  2.3, 

and  (e)  3.6  percent,  Gates  1977) 


CAVITATION  NUCLEI,  BUBBLE  MECHANICS  AND 
SCALING 

Nuclei  and  Their  Measurements 

Possibly,  no  other  subject  in  cavitation  inception  research  has  proved 
to  be  so  exacerbating  as  the  nucleation  sources  of  cavitation  or  "nuclei." 
The  many  theories  proposing  their  nature,  size,  and  occurrence  are  all  ex¬ 
tensively  reviewed  in  the  books  by  Knapp  et  al.  (1970)  and  Robertson  and 
Wislicenus  (1969).  Significant  advances  have  been  made  since  that  time  in 
the  detection  of  (rather  large)  nuclei  on  a  number  of  fronts;  these  are  re¬ 
viewed  first  by  Morgan  (1972),  and  then  later  by  Acosta  and  Parkin  (1975) 
and  do  not  need  to  be  further  elaborated  here  except  for  some  recent  find¬ 
ings  and  some  conceptual  notions.  The  basic  idea  is  rather  simple;  there 
must  be  some  practical  mechanism  for  the  phase  change  to  occur  in  engineer¬ 
ing  applications,  otherwise  the  large  tensions  predicted  for  pure  liquids 
far  removed  from  the  critical  point,  would  be  commonly  observed.  As  ex¬ 
plained  in  these  references,  nuclei  can  exist  as  small  gas-vapor  bubbles 
(perhaps  with  the  interface  contaminated  with  a  surfactant)  or  a  solid 
particle.  It  has  become  common  to  think  of  solid  particles  as  good  candi¬ 
dates  for  the  nuclei  because  of  the  theoretical  possibility  of  stabil¬ 
izing  a  gas  pocket  within  a  solid  particle.  Indeed,  Keller's  recent  work 
makes  a  case  for  these  "pore"  nuclei  in  his  own  facility  in  a  very  in¬ 
fluential  publication  (Keller  1973). 

The  measurement  techniques  have  been  both  optical  (single  point 
scattering,  see  Keller  1973)  and  holographic  photography,  (see  Peterson 
1972,  Peterson  et  al.  1975),  and  acoustic  and  indirect  via  the  Coulter 
counter.  The  acoustic  method  measures  only  bubbles;  both  the  Coulter 
counter  and  single  point  light  scattering  techniques  react  to  bubbles  and 
particulates.  As  in  photomicrography,  holographic  photography  is  limited 
in  resolution  but  the  entire  cross-section  of  tunnel  fluid  can  be  examined 
at  leisure  for  the  form  and  number  of  nuclei  greater  than  about  10  micro¬ 
meters  diameter.  (We  are  now  discussing,  of  course,  freestream  nuclei;  it 
is  also  possible  for  the  surface  of  the  body  itself  to  be  the  site  of 
nucleation,  just  as  it  is  in  most  boiling  situations.) 


From  various  sources,  we  summarize  (Figures  31  and  32)  the  number  of 
density  distributions  down  to  about  10  micrometers  diameter.  It  is  clear 
that  these  concentrations  can  vary  by  several  orders  of  magnitude;  several 
different  types  are  probably  represented  on  these  figures.  The  scattering 
data  contain,  undoubtedly,  both  solid  particles  and  microbubbles,  although 
for  the  larger  sizes  it  is  clear  from  the  reports  (Arndt  and  Keller  1976 
and  Keller  and  Weitendorf  1976)  that  air  bubbles  are  the  main  contributor. 
Gavrilov  (1970)  and  Medwin  (1977)  measure  microbubbles  acoustically.  Both 
particulates  and  bubbles  are  seen  in  the  holographic  studies  of  Peterson 
(1972),  Peterson  et  al.  (1975),  Gates  and  Bacon  (1977),  and  Katz  (1978). 

From  the  latter  we  present  an  example  showing  how  it  is  possible  to  readily 
detect  particu  lates  in  the  holographic  reconstruction  process  (Figure  33). 
We  could  add  to  these  plots  additional  data  from  the  St.  Anthony  Falls 
Hydraulic  Laboratory  (e.g.,  Schiebe  1969),  of  ocean  particulates  measured 
by  the  Coulter  counter  (Peterson  1974)  and  similar  measurements  by  the 
Naval  Ocean  Systems  Center.*  These  findings  all  fall  within  the  bands  of 
data  of  Figure  31,  and  they  all  have  about  the  same  functional  behavior 
for  the  small  sizes  (even  to  less  than  one  micrometer  for  the  particulates). 
This  apparent  uniformity  seems  quite  remarkable. 

Now  that  microbubble  and  particulate  nuclei  distributions  are  be¬ 
ginning  to  appear  more  frequently  in  the  hydrodynamic  literature,  one  may 
well  ask  just  how  these  data  are  to  be  used  quantitatively  in  cavitation 
inception  studies  and  in  flows  with  more  fully  developed  cavitation.  The 
several  orders  of  magnitude  difference  in  number  density  seen  in  Figure  31 
suggest  that  there  may  be  some  considerable  differences  to  be  expected  in 
cavitation.  We  have  in  fact  seen  such  differences  appearing  on  inception 
forms  already  in  Figure  6;  it  seems  certain  that  many  of  the  differences 
seen  there  are  due  to  the  nuclei  content.  It  is  also  clear  that  the 
Caltech  HSWT  has  a  very  low  bubble  nuclei  population  (Figure  32),  and  that 
the  numerous  particulates  of  that  facility  do  not  serve  as  nuclei.  Other 
investigators  have  shown  the  vital  importance  of  microbubble  supply  (but 
unquantified)  in  the  unsteady  flows  of  model  propeller  wakes  (e.g., 

Albrecht  and  Bjorheden  1975,  Noordzij  1976).  There  are  other  interactions 


*T.G.  Lang,  private  communication. 
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Figure  32  -  Nuclei  Distributions  Measured  by  Holography  in  the  Caltech 
LTWT  (All  Microbubbles)  and  HSWT  (Essentially  All  Solid  Particles) 

(Gates  and  Acosta  1978) 


between  bodies,  such  as  between  a  ship  hull  and  unsteady  propeller  cavi¬ 
tation,  where  the  nuclei  content  affects  the  gas  void  fraction  of  the 
medium  directly,  a  different  kind  of  problem  to  which  we  now,  momentarily, 
defer. 

To  make  full  use  of  the  nuclei  distribution  information,  we  need  an 
appropriate  physical  understanding  of  just  how  the  nuclei  themselves  lead 
to  cavitation  in  the  flows  of  interest.  (Direct  observation  of  this  process 
is  still  lacking.)  And  of  course,  one  would  like  to  predict  from  labora¬ 
tory  tests  (given  full  control  over  all  laboratory  variables)  what  cavita¬ 
tion  inception  indices  and  forms  might  be  on  previously  untried  prototype 
bodies  and  conditions.  We  think  this  is  still  not  quite  possible  except 
for  a  certain  kind  of  cavitation,  namely  travelling-bubble  cavitation.  To 
appreciate  the  mechanism  of  this  form  of  bubble  growth,  we  digress  briefly 
to  review  the  key  features  of  bubble  mechanics. 

Some  Bubble  Mechanics 

We  lean  heavily,  here,  on  the  monograph  of  Knapp  et  al.  We  are  con¬ 
cerned  with  the  growth  of  a  small  bubble  of  initial  radius,  Ro>  immersed  in 
a  pressure  field,  poo(t) ,  i.e.,  the  pressure  far  away  is  time  dependent. 

The  vapor  pressure,  p^,  depends  on  the  temperature  of  the  bubble  wall,  T(R), 
and  there  may  be  a  gas  partial  pressure,  p^.  We  need  the  equation  for  the 
motion  of  the  bubble  radius,  R(t),  which  is  (from  Knapp  et  al.  (1970)) 

+  !  *2  ■  5  [w4u  I  -  tt  -■>.]  <12> 

where  p  is  the  liquid  density,  y  the  viscosity,  and  ag  the  surface  tension 
and 

p  =  p  (T (R) )  +  p  (13) 

cav  v  g 

The  partial  pressure  of  the  gas  may  be  expressed  as 
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_  K(T) 
R3 


where  K,  proportional  to  the  mass  of  gas  within  the  bubble,  may  depend  upon 
the  temperature.  At  equilibrium,  the  bubble  radius,  R,  is  stationary,  say 
R  =  Rq.  Then  Equation  (12)  is 


o  =  -  (p„-pv>  -  jr  +  3 

o  R 

o 


Now  we  may  ask  whether  or  not  this  equilibrium  is  stable.  To  answer  this 
question  we  put,  in  the  usual  way. 


R  =  Ro  (l+e(t)) 


where  e(t)  is  supposed  to  be  a  small  number,  and  linearize  Equation  (12) 
to  get 


-6  +  ^£  +  jl  3K.ri 

R2  PR2Ir3  Ro 
o  o  \  o 


P°°  Pv  _  2°S  K 

P  _PRo  pR3 

o 


But,  the  right  hand  side  is  identically  zero.  We  have  now  the  equation  of 
a  single  degree  of  freedom  oscillator  and  we  expect  that 


e  =  £jfit 

o 


SI  being  the  frequency.  We  readily  find  that 


a  -  ± 

R 

O 


ks  -  nr  >  v  =  u/p 

R 

o 


These  conditions  may  be  termed  the  critical  ones  and  the  value  of  Rq,  so 
found  for  a  given  gas  content  of  bubble  (or  value  of  K)  the  critical  radius, 
R  .  Figure  34  shows  plots  of  these  equilibrium  radii  and  the  locus  of 
points  (Equation  (22))  separating  stable  and  unstable  regions.  Some  of  the 
numbers  from  Equations  (21)  and  (22)  are  interesting,  and  for  reference  we 
summarize  a  few  in  Table  2. 
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Figure  34  -  Equilibrium  Radius  for  Constant  Gas  Content  K  (Arbitrary  Units) 
Showing  Stable  and  Unstable  Regions  (Adapted  from  Knapp  et  al.  (1970)) 


TABLE  2  -  SOME  VALUES  OF  NATURAL  FREQUENCY  AND 
CRITICAL  RADIUS  FOR  BUBBLES 


R 

0 

8 

o 

(Poo-Pv  =  0-1  atm) 

(p  -p  ) 

K°°  Vent 

10“6  m 

2  x  106  Hz 

1.21  atm 

-5 

5 

10  m 

1  x  10  Hz 

0.12  atm 

-4 

A 

10  m 

1  X  10  Hz 

0.01  atm 

Before  continuing  on  to  consider  how  these  results  might  be  used  in 
interpreting  cavitation  inception,  we  pause  briefly  to  consider  the  vapor 
pressure  term,  p^,  and  how  it  may  be  expected  to  change  during  bubble 
growth.  The  basic  idea  is  that  bubble  growth  implies  evaporation  of  liquid 
and,  hence,  cooling  of  the  bubble  wall.  This  cooling  depresses  the  vapor 
pressure  and,  hence,  tends  to  reduce  the  pressure  difference  available  for 
growth.  Rough  ideas  of  the  magnitude  of  the  "thermal  effect"  are  easily 
obtained  following  physical  arguments  due  to  Plesset  (1949).  The  pressure 
difference  terms,  on  the  right  hand  side  of  Equation  (12),  are  split  up  as 
follows  (see  also  Holl  and  Kornhauser  1970) 

Pv<T)-Pco<t>  =  Pv(Too)-P00(t=0)  +  poo(t“0)-poo(t)  +  Pv(T)-pv(TJ  (23) 

where  (°°)  means  far  from  the  body,  and  T  is  the  temperature  of  the  liquid 
at  the  bubble  wall.  The  vapor  pressure  difference  is  estimated  from  the 
Clausius-Clapeyron  equation  to  be 

PV(T)  ~  PV(TJ  =  T“  Pv<T-TJ  <2*> 

cr 

where  is  the  density  of  the  vapor  and  £  is  the  latent  heat.  The  heat 
flow  supporting  the  evaporation  is  (approximately) 
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47TR2k  (T 

c  a 


-T) 


kc>  here,  being  the  conductivity,  and  6  a  heat  diffusion  distance.  This 
heat  flux  balances  the  evaporation  requirement 


C 


(26) 


and  Plesset  assumes  6  =  ✓’atT,  a  being  the  thermal  diffusivity.  Bringing 
these  together  we  have 


PV<T)-PV(T00) 

P 


FF^ 


The  other  pressure  differences  are  expressed  in  terms  of  the  conventional 

cavitation  number  and  pressure  coefficient.  To  express  the  scale  of  the 
*  * 

problem  let  t  =  t  L/U^,  where  t  is  a  dimensionless  time,  and  L  and  are 
a  reference  length  and  speed,  respectively.  We  may  then  write  the  right 
hand  side  of  Equation  (12)  (neglecting,  now,  the  viscous  term)  as 


RHS 
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(27) 


Two  interesting,  limiting  cases  are  immediately  evident:  first, 

neglecting  the  thermal  effect  entirely,  we  estimate  the  asymptotic  growth 

rate  for  an  imposed  constant  value  of  (a  +  c  ) .  Then  for  R  >>  R  we  have 

p  o 

the  limiting  "inertial"  growth  rate 

1 

”  Uo  (-T  <c,  +  v)2  <28) 
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a  constant.  Second,  we  may  now  imagine  that  the  thermal  effect  is  large, 
so  that  the  bubble  growth  is  slow.  Then,  the  inertial  terms  on  the  left 
hand  side  of  Equation  (12)  may  be  imagined  to  be  small.  Again,  for  R  »  R 

i 

bubble  growth  is  then  in  thermal  equilibrium  with  a  rate 


*th  =  -  2- 


(o  +  cp) 


(29) 


These  rates  can  be  quite  different.  In  the  first  case  we  say  that  we 
have  "cavitation,"  and  in  the  second,  "boiling"  is  said  to  occur.  The 
ratio  of  the  two  rates  may  be  expressed  as 


th 


_  2 

V3 


/  t*  /  M  c2  /LUo  .  i_ 

y  “(a  +  cp)  \p  /  cpToo  y  a  y2 


(30) 


and  for  boiling  to  occur,  the  thermally  controlled  rate  must  be  much  less 
than  the  inertial  one,  i.e., 


»  1 


th 


(31) 


This  has  led  Brennen  (1973)  to  identify  the  dimensional  group  of  terms 


as  the  "thermal  parameter"  of  importance  and  thence,  if 


(32) 
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z»y: 


-(o  +  y 


(33) 


boiling  occurs.  For  many  applications  in  naval  hydrodynamics,  we  may 
imagine  that  the  term  -(a  +  c  )/t  *  1  so  that  Brennen's  \  function  can 
be  made  dimensionless  with  the  result  that  if  the  term 


(34) 


boiling  is  said  to  occur.  In  this  group  of  terms,  E^  is  the  Eckert  number, 
and  P£  is  the  usual  Peclet  number,  both  familiar  in  heat  transfer.  An 
example,  here,  may  illustrate  the  point.  We  consider  water  with  reference 
speed  of  10  m/sec  and  reference  length  of  3  cm  to  form  Table  3. 


TABLE  3  -  COMPARISON  OF  THERMAL  PARAMETERS  FOR  WATER 

U  =  10  m/s,  L  =  3  cm 
o 


T  =  20°  C 

260°  C 

360°  C 

Et  =  1.4  x  10-9 

-4 

7  x  10’ 

5  x  io-3 

E.  /ST  =  1.5  X  10-6 

k  e 

1.5  x  io-6 

1.5  x  io'6 

— ►  Cavitation 

Boiling 

Boiling 

We  see  then,  that  these  thermal  effects  are  not  important  for  hydro¬ 
dynamics,  but  may  be  for  power  machinery.  Brennen  (1973)  includes 
comparisons  with  many  other  fluids,  and  concludes,  interestingly,  that 
liquid  hydrogen  almost  always  boils  rather  than  cavitates. 
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Scaling  Ideas 

With  these  results  from  bubble  mechanics  in  hand  the  main  features  of 
travelling-bubble- inception  scaling  scenarios  become  apparent:  A  micro¬ 
bubble,  or  nucleus  of  some  size,  may  be  imagined  to  be  swept  through  the  low 
pressure  region  of  a  flow.  The  small  nuclei  have  a  very  high  natural  fre¬ 
quency  and  so  react  quasi-statically  to  the  imposed  pressure  change.  It 
may  happen  that  the  critical  pressure  difference  is  exceeded  for  a  certain 
freestream  microbubble;  then  it  follows  that  there  will  be  a  period  of 
exponential  growth  which,  if  of  sufficient  duration,  will  ensure  a  visible 
macroscopic  bubble.  Indeed  many  workers  (e.g.,  Johnson  and  Hsieh  1966, 
Schiebe  1972,  Silberman  1975)  are  concerned  only  with  that  portion  of  the 
fluid  flow  at  below  the  critical  pressure,  since  the  subsequent  bubble 
growth  is  so  rapicj.  There  is  an  additional  possibility  that  turbulent 
pressure  fluctuations  in  transition  locations  may  actually  excite  nuclei 
there  to  unstable  growth  by  resonance.  In  any  case,  from  such  arguments  it 
is  evident  that  the  number  of  travelling-bubble  cavitation  "events"  will  be 
proportional  to  the  number  of  microbubbles  that  can  grow  exponentially  or 
explosively;  this  then,  depends  upon  the  nuclei  density  distribution,  the 
pressure  difference  (p^-p  ) ,  and  the  particular  body. 

Inception  for  travelling-bubble  cavitation,  then,  is  determined  by 
the  rate  of  such  events  for  a  particular  body.  This  concept  of  "event 
counting"  as  previously  mentioned  has  been  followed  vigorously  by  the  St. 
Anthony  Falls  Hydraulic  Laboratory  (e.g.,  Silberman  1975,  Schiebe  1972). 

In  fact,  these  authors  turn  the  process  around  to  infer  the  number  density 
distribution,  N(R),  by  the  number  of  such  events  on  "standard"  bodies. 

Apart  from  this,  their  measurements  of  events  show  impressive  correlation 
with  their  acoustically  determined  nuclei  distributions.  These  facilities 
contain  many  nuclei  at  the  100  micrometer  level  so  that  only  a  slight 
pressure  difference  (poo~Pv)  is  needed  to  establish  criticality  (see,  e.g., 
Table  2).  But  travelling-bubble  cavitation  is  rare  in  the  Caltech  HSWT, 
and  tensions  of  up  to  1/2  atmosphere  are  readily  sustained  in  flows  past 
bodies  at  minimum  pressure  points.  From  this  we  conclude  that  freestream 
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nuclei,  larger  than  about  5  micrometers  diameter,  are  not  very  abundant  in 
this  facility,  a  conclusion  supported  by  the  data  of  Figure  32.  Further¬ 
more,  it  points  out  that  the  range  of  sizes  important  in  cavitation  testing 
work  should  extend  downwards  at  least  that  far. 

There  is,  then,  a  direct  influence  of  the  nuclei  density  distribution 

on  the  rate  of  travelling-bubble  cavitation  events;  when  there  are  very  few 

freestream  nuclei,  these  events  will  of  course  be  correspondingly  few. 

And,  as  we  have  seen,  appreciable  tensions  can  be  readily  achieved,  even  in 

comparatively  "dirty"  tunnel  water.  From  the  physical  evidence  of  tests 

in  several  tunnels,  attached  forms  of  cavitation  may  occur  at  inception 

even  in  the  absence  of  laminar  separation.  Again,  we  return  to  the  Caltech 

facilities  in  Figure  35  to  show  these  different  types  of  inception  on 

identical  test  bodies.  Not  surprisingly,  it  was  found  that  the  inception 

indices  differ  markedly  with  travelling-bubble  cavitation,  showing  the 

higher  values  (although  still  considerably  below  |-c  |).  In  fact,  in 

Pmin 

the  HSWT  there  was  a  closer  correlation  with  -c  than  any  other  parameter 

Ptr 

(Gates  and  Acosta  ibid).  That  attached  forms  of  cavitation  can  sometimes 
occur  first  even  on  nonseparating  bodies  shows  how  difficult  it  is  to  have 
a  single  universal  threshold  measure  of  inception. 

We  can  now  recognize  that  a  single  kind  of  phenomenon  that  can  lead 
to  a  uniform  scaling  law  in  cavitation  inception  does  not  occur.  It  is 
worthwhile,  however,  to  draw  together,  in  Table  4,  the  following  "expecta¬ 
tions"  or  trends  in  inception  as  we  now  understand  it  based  on  water 
tunnel  tests  (from  Arakeri  and  Acosta  1979).  It  would  be  extremely 
valuable  to  extend  these  observations  and  even  trends  to  prototype 
Reynolds  numbers;  this,  however,  must  await  a  future  day. 

FULLY  DEVELOPED  CAVITATION  IN  INTERNAL 
FLOWS 

We  now  leave  inception  to  consider  more  developed  states  of  cavita¬ 
tion.  This  is  a  subject  for  which  vast  literature  exists.  Our  present 
purpose  is,  however,  a  more  modest  one;  namely  to  consider  developed 
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SCHIEBE  BODY,  LTWT.  24  ft. /sec. 
D  =  2  inches,  a  =  0.52 
POLYOX  INJECTION 


SCHIEBE  BODY,  HSWT,  46  ft./sec. 

D  =  2.0  inches,  o  =  0.41 ,  S/D  =  0.47 

Figure  35  -  Cavitation  on  the  Same  (Schiebe)  Body  in  the  Caltech  LTWT  and 
the  Caltech  HSWT  (Gates  and  Acosta  1978) 


TABLE  4  -  TYPES  OF  CAVITATION  INCEPTION  PHENOMENA 


Nuclei  Supply 

Type  of  Flow 
on  Body 

Expected  Inception 

Forms 

COPIOUS 

No  Separation 

Travelling  bubble  from  c 
(e.g.,  Silberman  1975).  Pmin 

Event  counting. 
a  <  -c 

Pmin 

COPIOUS 

Laminar  Separation 

Mixed  travelling-bubble  and 
band  cavity  (sometimes  fore¬ 
stalled  by  myriads  of  bubbles) 

FEW 

No  Separation 

Attached  cavities, 
a  -  -c 

Ptr 

(Arakeri  et  al.),  wedges,  fingers, 
patches  often  seen  up  to  position 
of  c 

p  . 
min 

INTERMEDIATE 

No  Separation 

Travelling-bubble  cavitation/ 
attached  forms, 

(e.g.,  Huang  and  Santelli  1977) 

o  =  -  c 

Ptr 

INTERMEDIATE 

AND 

FEW 

Laminar  Separation 

Band  and  sheet  cavities.  Glassy 
leading  edges 

a  -  -  ( c  +Ac  ( t ) ) 

Ps  P 

(e.g.,  Huang  and  Peterson  1976). 

NOTE:  Desinent  cavitation  and 

gaseous  inception  are  not  included. 

cavitation  in  internal  flows.  These  are  the  flows  that  take  place  in 
water  tunnels,  fluid  machinery,  and  the  duct  flows  of  various  types  that 
arise  in  large  hydraulic  systems  and  the  propulsion  systems  of  ships  (see, 
e.g.,  Narita  and  Kunitake  1977,  Yamaguchi  1977).  As  in  cavitation  in¬ 
ception,  different  forms  of  developed  cavitation  occur,  and  for  internal 
flows  it  is  convenient  to  lump  them  into  three  types,  as  in  Table  5: 


TABLE  5  -  TYPES  OF  INTERNAL  FLOW  CAVITATION 


1. 

Attached 

two-  and  three-dimensional 

surface 

cavities. 

2. 

Tip-vortex  and  tip-clearance  flow 
cavitation. 

3. 

"Bubbly" 

cavitation  and  two-phase  flow. 

Although  all  of  these  forms  are  commonly  seen  in  water  tunnel  tests,  type  2 
is  most  often  associated  with  flow  in  fluid  machines,  such  as  a  propeller 
or  pump.  Often  the  tip-clearance  flow  in  a  pump  is  sufficiently  strong  to 
cause  a  significant  "backflow"  resulting,  thereby,  in  a  very  complex  inflow 
to  the  pump. 

Attached  cavity  flows  have,  of  course,  received  much  attention  because 
of  the  powerful  analytic  tools  that  can  be  brought  to  bear,  as  well  as  the 
intrinsic  importance  of  these  flows  (see,  e.g.,  Wu  1972).  These  same 
techniques  may  be  used  to  explain  tunnel-caused  effects  in  steady  flow 
which  may  be  termed  "wall  effects"  (see,  e.g.,  the  review  by  Baker  1977). 
When  the  void  fraction  of  liquid-gas-vapor  flows  is  appreciable  it  is 
customary  to  call  these  "two-phase  flows"  (see,  e.g.,  the  text  by  Wallis). 
Indeed,  much  effort  has  been,  and  is  being,  expended  on  this  subject  as  it 
is  a  key  issue  in  the  understanding  of  the  thermo-hydraulics  processes  in 
nuclear  power  plants.  But  many  cavitating  flows  in  naval  hydrodynamics 
consist  of  dispersed  cavitation  bubbles  in  which  the  void  fraction  is  well 
below  one  percent  and  these  flows  we  may  term  bubbly  cavitation.  This 
kind  of  cavitation  is  common,  even  dominant,  in  cavitating  pumps  and,  as  we 
have  seen,  is  also  a  highly  visible  feature  of  some  water-tunnel  flows. 
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A  special  feature  of  internal  flows  with  cavitation  is  the  coupling 
of  the  fluid  motion  in  various  parts  of  a  system  by  the  displacement  of  a 
changing  volume  of  cavitation.  This  unsteadiness  may  arise  deliberately 
as  a  part  of  a  test  procedure  (e.g.,  a  study  of  propeller-wake  inter¬ 
action),  or  it  may  happen  that  the  cavitation  itself  is  inherently  unsteady. 
(Partial  cavitation  of  about  three-quarters  chord  length  provides  such  an 
example.)  Because  of  the  coupling  of  the  flow  with  unsteady  cavitation, 
there  is  an  interaction  between  the  internal  flow  circuit  and  the  basic 
flow  past  the  cavitating  body  itself.  This  interaction  may,  in  fact,  so 
control  the  whole  process  that  the  original  basic  flow  may  be  entirely 
masked  (as  for  example  in  an  unsteady  hydrofoil  test).  The  confining 
effect  of  a  tunnel  has  been  in  fact  a  serious  deterrent  to  measurement  of 
unsteady  forces  on  cavitating  hydrofoils  in  the  past.  We  need  ways,  then, 
to  correct  for  the  "unsteady"  internal  effect,  as  well  as  for  the  better- 
understood  steady  wall  effect. 

It  may  happen  that,  because  of  cavitation  and  circuit-cavitation 
interactions,  an  inherently  steady  or  stable  cavitation  flow  may  exhibit 
self-excited  oscillations.  Such  oscillations  occur  frequently  in  pump 
applications  or  in  more  complex  systems,  such  as  pump-fed  liquid  fueled 
rocket.  This  is  the  basis  of  the  notorious  "Pogo"  longitudinal  oscilla¬ 
tions  of  large  liquid  booster  rockets.  In  fact,  some  of  these  Pogo 
oscillations  have  reached  acceleration  levels  of  50  g's*  and  it  is 
customary  for  all  liquid  fueled  rockets  to  suffer  some  extent  of  this 
instability.  Cavitating  pumps  are  usually  thought  of  primarily  in 
connection  with  space  applications.  However,  the  types  of  designs  and 
performance  demands  of  modern  high-speed  naval  craft,  such  as  the  PHM 
hydrofoil  and  surface  effect  ship  propulsion  systems,  are  very  similar  to 
those  for  large  booster  propulsion.  But  it  is  not  necessary  to  go  so  far 
afield  to  find  situations  where  unsteady  cavitation-system  coupling  presents 
a  problem  of  interpretation.  A  good  example  is  afforded  by  the  unsteady 
hull  pressure-propeller  cavitation  measurements  of  Keller  and  Weitendorf 
(1970).  In  their  tests,  the  pressure  at  several  points  on  a  wall  adjacent 
to  a  cavitating  propeller  was  measured  in  a  closed  tunnel  as  a  function  of 


*Dr.  S.  Rubin,  Aerospace  Corp. 


72 


undissolved  void  fraction  (10  ^  to  10  ^).  The  pressure  was  found  to  vary 

strongly  with  both  frequency  and  void  fraction  when  operating  in  a  non- 

-4 

uniform  wake  except  at  very  high  void  fractions  (10  )  where  the  results 

appeared  to  be  independent  of  frequency.  These  interesting  observations 
show,  again,  the  importance  of  the  presence  of  microbubbles  on  the  extent 
of  cavitation,  and  they  raise  the  question  of  tunnel  interaction  through 
the  confining  effect  of  the  surrounding  walls  on  a  dynamic  experiment. 

We  turn,  now,  to  the  immediate  objectives  of  the  present  section.  We 
first  touch  upon  the  available  analytic  methods  that  can  be  used  to  esti¬ 
mate  these  interaction  effects  for  the  types  of  cavitation  mentioned, 
and  we  then  discuss  some  recent  experimental  work  aimed  at  characterizing 
unsteady  flows  through  cavitating  pumps. 

Analytic  Methods  for  Internal  Cavitation 

We  refer  now,  briefly,  to  Table  5  for  the  types  of  cavitation  to  be 
considered.  We  are  concerned  here,  primarily,  with  unsteady  problems.  It 
is  clear  from  previous  experience  that  attached  unsteady  cavities  can  be 
dealt  with  on  the  basis  of  linearized  potential  flow  (Wu  1972).  There 
are,  however,  still  difficult  problems  for  the  proper  formulation  of 
unsteady,  growing,  attached  cavities  with  realistic  geometries.  The  tip 
vortex  flow  is  also  complex;  dynamic  treatment  of  these  flows  appear  to 
be  lacking.  Pumps  and  propellers  reveal  attached  cavities,  as  well  as 
bubbly  tip-clearance  flows  and  inflows.  The  global  properties  of  a 
volume  of  such  bubble  flow  can  be  established  with  bubble  mechanics  for 
unsteady  pressure  fluctuations.  But  the  concentration  of  bubbles  in  such 
flows  must  be  small,  otherwise  the  full  two-phase  fluid  dynamic  equations 
for  the  mixture  of  cavitating  bubbles  and  liquid  must  be  used.  These 
equations  are  evidently  not  yet  well  developed  analytically.* 

In  the  following  sections,  we  sketch  a  sample  problem  intended  to 
illustrate  some  of  these  dynamic  features  for  a  pump. 


*Private  communication  with  Prof.  L.  van  Wijngaarden. 
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The  Coupling  Problem  for  a  Pump 

We  sketch  the  blades  of  an  axial  inducer  together  with  the  types  of 
cavitation  in  Figure  36.  We  imagine  that  the  upstream  pressure  and  inlet 
velocity  are  fluctuating.  We  anticipate  that  the  outlet  quantities  of  the 
pump  will  similarly  fluctuate.  There  is  upstream  structure  and  downstream 
structure  to  input  and  receive  the  unsteady  flow  through  the  pump.  Our 
coupling  problem  here  is  to  provide  the  effect  due  to  the  insertion  of  the 
pump  or  other  cavitating  device  on  these  up  and  downstream  flows.  It  is 
natural,  then,  to  think  of  the  pump  as  having  a  "transfer  function"  whereby 
periodic  input  disturbances  are  connected  to  output  ones.  An  extended 
system  can  then  be  treated  by  connecting  its  various  transfer  functions 
appropriately.  In  what  follows,  some  components  of  such  a  transfer 
function  are  estimated  from  the  standard  reference  problem  in  turbo¬ 
machinery  analysis,  the  flow  through  a  cavitating  cascade.  The  flow  model 
is  that  of  an  attached  blade  cavity  with  a  fluctuating  axial  inlet  flow 
velocity.  The  cavity  is  assumed  to  be  at  constant  pressure;  its  boundaries 
and  volume  change  in  some  relationship  to  the  inlet  flow  perturbations. 

To  focus  on  the  effects  of  cavitation,  the  wetted  part  of  the  cascade  is 
taken  to  be  infinitely  long.  This  representation,  at  first  sight,  highly 
artificial,  is  actually  a  good  one  for  the  partial  cavitation  within  the 
inducer  portions  of  cavitating  pumps.  A  sketch  of  the  cascade  is  given  in 
Figure  37a  and  linearized  boundary  conditions  in  Figure  37b.  This 
linearized  boundary  value  problem  implicit  in  Figure  37b  is  capable  of 
formal  solution  in  an  auxiliary  plane  where  the  entire  blade  and  cavity 
geometry  appears  on  the  real  axis  by  the  methods  outlined  in  Wu's  review. 

We  need  not  concern  ourselves  herein  with  these  details  except  for  several 
basic  features  of  the  solution.  There  must  be  a  fluctuating  pressure  gra¬ 
dient  approaching  the  cascade  because  of  the  fluctuating  inflow  velocity. 
There  will  be  a  similar  fluctuating  downstream  pressure  gradient,  possibly 
with  a  different  amplitude  and  phase,  because  of  the  changing  cavity 
volume.  We  can  anticipate  that  the  actual  pressure  will  be  of  the  form 
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Figure  37a  -  Physical  Plane 


Figure  37b  -  Mathematical  Plane  (with  linearized  boundary  conditions) 


Figure  37  -  Representation  of  Unsteady  Cavitating  Flow  Through  a  Cascade 
of  Inducer  Blades  Together  with  Linearized  Boundary  Conditions 

(Kim  and  Acosta  1975) 
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downstream.  Here  w  is  the  frequency,  j,  the  imaginary  unit  of  time, 
u^  and  u2>  the  fluctuating  velocities  in  Figure  37b,  and  x,  the  distance 
along  the  chord.  The  first  term  may  be  seen  to  be  the  pressure  gradient 
due  to  the  "mass  oscillations"  of  the  up  and  downstream  flows.  The  terms 
p^s  and  p2s  are  the  steady  pressures  across  the  cascade  in  the  absence  of 
motion,  and  p^  and  p>2  are  the  additional  complex  fluctuating  pressures 
caused  by  the  presence  of  the  cascade. 

These  additional  pressures  are  what  are  needed  to  account  for  the 
insertion  of  the  cascade  between  the  up  and  downstream  flow  regions.  It 
turns  out  that  these  additional  pressures  can  be  simply  related  to  the 
unsteady  axial  fluctuations  by  a  linear  equation  of  the  form 


where  z  is  a  2  x  2  matrix  with  complex  coefficients.  This  matrix,  we  can 
now  see,  is  the  transfer  function  for  the  cascade.  Before  discussing  some 
of  the  terms  of  this  transfer  function,  we  need  to  mention  some  of  the 
difficulties  in  finding  practical  evaluations  of  the  implied  formal 
solution. 

Even  in  the  linearized  free  streamline  theory,  we  need  to  be  concerned 
about  the  physical  significance  of  the  cavity  model.  It  has  become 
acceptable  to  think  of  steady  flow  cavity  models  as  having  an  ultimate 
"wake"  to  account  for  the  momentum  deficiency  of  the  cavitating  drag. 
Unsteady  cavitation  models  have  not,  however,  received  as  much  attention. 
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In  water  tunnels,  we  observe  that  attached  cavities  grow  and  shrink  in  re¬ 
sponse  to  system  pressure  changes.  To  the  eye,  at  least  for  slow  changes, 
the  cavity  appears  to  be  a  closed  body  of  changing  shape  and  length.  It 
is  actually  a  material  surface.  In  the  linear  theory  the  coordinate,  Y, 
of  the  end  of  the  cavity,  is 


Y(t)  = 


(38) 


where  SL (t )  is  the  fluctuating  end  of  the  cavity,  v  is  the  perturbation 
velocity  component  normal  to  the  cavity,  and  is  the  reference  speed. 

For  steady  flows,  Y  =  0  is  Tulin's  closed  cavity  model;  Y  ^  0  corresponds 
to  an  open  wake  model.  For  unsteady  flows,  Equation  (38)  is  actually  quite 
a  difficult  integral  to  evaluate,  even  in  simple  flows.  Most  cf  these 
problems  were  foreshadowed  by  Parkin's  work  on  unsteady  hydrofoils 
(Parkin  1957).  As  it  has  turned  out,  most  work  in  the  field  since  then  has 
avoided  direct  use  of  Equation  (38),  and,  because  of  these  formal  diffi¬ 
culties,  other  approaches  have  been  adopted.  Among  these,  two  are  well- 
known;*  Guerst  permitted  no  change  in  cavity  volume,  whereas  Leehey  assumed 
that  the  end  of  the  cavity  remained  at  its  steady  value  (thereby  circum¬ 
venting  Equation  (38)),  but  allowed  the  cavity  volume  to  change.  This  ap¬ 
proach  was  also  used  later  by  Kim  and  Acosta  (1975)  to  evaluate  the  terms 
of  the  transfer  matrix,  z.  Still,  this  seemed  not  entirely  a  satisfactory 
resolution  because  it  was  clear  that  very  slow  oscillations  should  approach 
known  "quasi-steady"  values,  wherein,  the  actual  changes  in  cavity  length 
that  are  seen,  are  properly  modeled.  It  has  been  shown  more  recently 
(Acosta  and  Furuya)  that  application  of  Equation  (38),  with  Y  =  0  for  small 
frequency,  does,  in  fact,  lead  to  the  quasi-steady  limit  of  Tulin’s  closed 
cavity  model.  It  was  also  shown  that  the  f ixed- terminus  cavity  model  for 
unsteady  flow  had  a  rather  unnatural  singular  behavior  of  one  of  the 
coefficients  of  z.  To  appreciate  this  we  now  need  to  comment  on  each  of 
the  terms  of  the  cascade  transfer  function  z. 


*Refer  to  the  review  by  Wu  for  details. 


The  Transfer  Function 


The  sample  problem  of  this  section  is  in  the  form  of  the  flow  across  a 
moving  cascade.  The  static  pressure  rise  across  such  a  cavitating  cascade 
is  a  function  of  the  upstream  pressure  and  the  upstream  flow  velocity.  In 
the  absence  of  fluid  losses,  an  infinitely  long  cascade  has  no  change  in 
pressure  rise  with  inlet  pressure  for  a  steady  flow.  An  unsteady  flow  is  a 
different  matter,  however,  and  the  component,  z.^,  may  not  be  expected  to 
be  zero.  The  fluid  losses  in  steady  flow  are  due  primarily  to  mixing 
caused  by  the  drag  due  to  cavitation.  This  "momentum"  loss  increases  with 
decreasing  pressure,  and  thus,  we  may  expect  quasi-steady  cascade  experi¬ 
mental  data  to  show  a  positive  value  of  z^.  Pump  cascades,  as  sketched  in 
Figures  36  and  37,  create  a  static  pressure  rise  through  diffusion  of  the 
relative  velocity.  This  pressure  rise  decreases  with  increased  axial 
speed  (and  thereby  reduced  angle  of  attack)  so  that  on  a  quasi-steady  basis 
we  may  expect  the  term  z^  to  be  negative,  i.e.,  a  fluctuating  increase  in 
flow  should  result  in  a  fluctuating  decrease  in  pressure  rise.  Additional 
complex  terms  may  be  expected  in  unsteady  flow. 

We  have  repeatedly  focused  attention  on  the  volume  change  of  cavita¬ 
tion  with  pressure;  the  rate  of  volume  change  is  seen  to  be  proportional  to 
the  differences  of  down  and  upstream  speed  across  the  cascade.  From 
experience  with  steady  cavitating  flows,  we  expect  the  volume  of  cavita¬ 
tion  to  increase  with  a  decrease  in  upstream  pressure,  and  to  increase  with 
an  increase  of  angle  of  attack  or  a  decrease  in  inlet  flow  velocity.  We 
expect  that,  for  vanishing  frequency,  the  terms  z ^  and  z ^  will  both  be 
negative  and  imaginary,  and  proportional  to  frequency,  as  this  is  the 
appropriate  behavior  in  the  quasi-steady  limit.  The  term  is  of 
particular  interest;  it  is  related  to  what  has  been  termed  the  "compliance" 
by  hydraulic  system  analysts.  The  compliance  is  the  negative  rate  of 
change  of  cavity  volume  with  pressure.  In  fact,  in  dimensional  units, 
is  equal  to  the  negative  of  this  compliance  times  jio  divided  by  the  pump 
inlet  area.  It  is  an  important  system  parameter,  since  the  period  of  an 
inlet  pipe  line  oscillation  is  proportional  to  the  square  root  of  the 


compliance  times  the  inertance  of  the  inlet  line.  It  is  the  compliance 
that  becomes  singular  for  the  fixed-terminus  cavity  model,  thereby  showing 
this  model  to  be  unsuitable  for  an  unsteady  internal  flow. 

Thus  far,  in  discussing  unsteady  flow  interactions,  we  have  focused  on 
the  cascade  idea.  This  was,  primarily,  because  the  natural  form  of  the 
solution  led  immediately  into  the  transfer  function  description,  more 
common  in  system  analysis.  As  mentioned  in  Table  5,  not  all  flows  of 
interest  are  attached  cavity  flows;  in  fact,  these  may  not  be,  dynamically, 
the  most  important  forms.  Observations  of  some  pump  inflows  reveal  a 
significant  number  of  microbubbles.  These  undergo  growth  and  collapse 
histories  within  the  impeller  which  are  modified  by  the  fluctuating  inlet 
pressure.  The  "compliance"  of  a  stream  of  microbubbles  undergoing  such  a 
growth  and  collapse  process  can  be  determined  from  bubble  mechanics.  This 
is,  in  fact,  what  Brennen,  in  a  most  influential  paper,  has  carried  out  in 
great  detail  (Brennen  1973),  using,  for  want  of  a  better  choice,  microbubble 
distributions  derived  from  water  tunnel  measurements.  The  results  are  most 
interesting  in  that  the  response  of  individual  bubbles  shows  that  the 
compliance  slowly  decreases  from  its  quasi-static  value  as  frequency  in¬ 
creases  with  a  gradually  increasing  phase  lag.  Thus,  in  the  present 
notation,  becomes  complex  with  increase  in  frequency.  We  defer  these 
interesting  results  for  comparison  later  with  cavitating  pump  applications. 

Application  to  Pumps 

Pump  engineers,  historically,  have  carried  out  two  types  of  steady 
tests;  namely,  a  performance  test  in  which  the  (total)  pressure  difference 
at  fixed  rotative  speed  is  measured  as  a  function  of  flow  rate,  and  then, 
a  cavitation  test  in  which  the  pressure  rise  at  fixed  speed  and  flow  rate 
is  measured  as  a  function  of  inlet  pressure.  It  is  customary  to  normalize 
the  total  pressure  rise  by  the  quantity  (q),  and  the  inlet  pressure  by  (q), 
where  q  is  the  dynamic  pressure  based  on  tip  speed.  The  axial  inlet 
velocity  is  normalized  by  the  tip  speed  (different,  thereby,  by  tt  from  the 
advance  ratio  of  propellers).  Curves,  representative  of  the  performance  of 
real  pumps,  are  sketched  in  Figure  38.  From  our  previous  discussion  of 
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unsteady  flows  across  cavitating  pump  cascades,  we  can  see  that  a  plausible 
form  of  the  transfer  function  for  a  pump  would  be  (with  reference  to  the 
notation  of  Figure  38) 


(39) 


the  circumflexes,  again,  refer  to  the  unsteady  fluctuations.  In  the  limit 
of  zero  frequency,  the  steady  pump  performance  curve  must  have  the  result 
that 


z1 1  -+  tan  a 

(40) 

z12  ->  tan  g 

for  infinitesimal  excursions  about  the  operating  point.  When  the  frequency 
is  not  zero  we  may  expect  that  z^  and  z^2  become  complex  and  that  z^  and 
z22  will  become  important 

The  steady  state  pump  test  then  provides  limiting  estimates  of 
and  z-^2*  Much  of  the  early  work  on  hydraulic  system  dynamics  involving 
cavitating  pumps  has  been  based  upon  the  idea  that  the  term  z^  is  pre¬ 
cisely  a  "compliance"  effect  and  that  z22  was  zero.  The  compliance  was 
then  determined  by  making  a  best  fit  to  field  data  of  observed  system 
oscillations  (see,  e.g.,  Rubin  1966).  It  was,  of  course,  not  necessary  to 
assume  z22  to  be  zero  except  for  simplicity.  Later,  Brennen  et  al.,  (1976) 
showed  from  limiting  quasi-steady  arguments  based  on  field  data,  that  the 
term  corresponding  to  z22  was,  indeed,  important,  and  that  this  term,  as 
well  as  could  be  estimated  from  quasi-steady  cavitating  cascade  theory. 

As  it  turned  out,  these  estimates,  even  when  the  effects  of  blade  thick¬ 
ness  were  included,  underestimated  the  results  of  field  observations  by  a 
fair  amount,  particularly  so  for  the  compliance  term.  It  had  been  shown 
previously  by  Brennen  (Brennen  1973)  that  the  compliance  effect  due  to 
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travelling  microbubbles  alone  could  account  for  much  of  the  observed  be¬ 
havior.  We  reproduce  here  in  Figure  39  his  summary  results  using  the 
water-tunnel  derived  estimates  of  microbubble  populations.  Plainly  al¬ 
though  much  remains  to  be  learned  about  the  cavitating  environment  of  these 
flows,  bubbly  cavitation  is  evidently  an  important  factor  in  them.  From 
all  of  this  work,  it  was  clear  that  new  experimental  methods  had  to  be 
devised  to  measure  the  transfer  function  of  a  pump. 

Experiments  in  Unsteady  Pump  Dynamics 

There  has  been  growing  concern  by  many  workers  with  the  proper 
representation  of  unsteady  pump  characteristics.  We  may  mention,  for 
example,  the  work  of  Kolesnikov  and  Kinelev  (1973),  Fanelli  (1972),  and 
Black  and  Santon  (1975),  who  are  all  concerned  in  various  respects  with 
this  problem.  Experimental  work  is  comparatively  rare,  however.  Among  the 
first  of  these  is  that  of  Ohashi  (1968)  and  later  Anderson  et  al.  (1971), 
who  carried  out  experiments  on  oscillating  flow  in  fully  wetted  pumps. 
There,  in  principle,  only  the  term  z ^  needs  to  be  measured;  this,  of 
course,  greatly  simplified  measurement  demands.  Rather  similar  problems 
have  encroached  into  related  fields  of  naval  hydrodynamics,  namely  the  re¬ 
sponse  of  lift  fans  of  surface  effect  vehicles  to  periodic  disturbances. 
There  too,  a  transfer  function  approach  is  taken  (Durkin  and  Luehr  1978). 
Durkin  and  uuehr  did  not  carry  out  measurements  of  the  transfer  function 
itself,  but  based  their  conclusions,  as  did  earlier  workers,  on  lumped 
parameter  system  model  tests.  It  would  appear  that  the  first  direct 
attempt  to  measure  the  transfer  function  of  a  cavitating  pump  system  is 
due  to  Ng  (1976).  The  problem,  as  described  there,  is  to  determine  the 
eight  unknown  coefficients  of  the  transfer  matrix  [z]  of  Equation  (39)  from 
measurements  of  fluctuating  inlet  and  outlet  pressures  and  volumetric  flow 
rates.  Any  given  hydraulic  system  with  a  fixed  method  of  excitation  will 
produce  a  unique  relation  between  input  and  output  quantities,  thereby 
providing  insufficient  data  to  determine  [z].  The  solution  devised  by  Ng 
was  to  provide  two  sources  of  excitation  to  a  test  loop,  these  sources 
being  locked  in  phase,  but  separated  by  an  isolation  section.  Then,  with 
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Figure  39  -  Cavitation  Compliance  K_  Estimated  from  Microbubble  Populations 
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Measured  in  Water  Tunnels  versus  Cavitation  Index  (Brennen  1973) 


fixed  operation  parameters,  variation  of  the  phase  between  the  sources 
provides  the  additional  sets  of  relations  needed  to  find  [z].  In  practice 
combinations  of  phase  and  amplitude  were  used.  The  excitation  consisted 
of  a  "siren"  like  valve  which  provided  a  variable  hydraulic  resistance  to 
the  flow  circuit. 

Another  notable  feature  of  the  Ng  experiments  was  the  use  of  laser 
velocimetry  to  measure  up  and  downstream  flows.  This  necessitated  certain 
flow  smoothing  and  conditioning  not  common  in  hydraulic  circuits.  These 
experimental  features  are  shown  in  Figure  40,  and  an  example  of  some  of  the 
first  experimental  results  are  shown  in  Figure  41  (Ng  et  al.  1976). 

This  rather  involved  figure  shows  two  sets  of  data,  one  basically 
noncavitating  and  the  other  cavitating.  The  noncavitating  results  show 
that  the  terms  z^,  z2\'  anc*  z22  are  sma^»  in  t^le  absence  of  cavitation 
and  elasticity  of  the  fluid  and  supporting  structure,  and  experimental 
error,  these  terms  should  be  strictly  zero.  The  term  z ^  is  seen  to  have 
a  negative  real  part.  This  corresponds  fairly  well  to  the  slope  of  the 
total  head  performance  curve,  the  angle  B  of  Figure  38a,  but  shows  some 
change  with  frequency.  The  term  z ^  also  has  an  imaginary  part  which  is 
equivalent  (within  experimental  accuracy)  to  the  inertia  of  a  certain 
length  of  fluid.  The  fully  wetted  data  are  rather  as  expected. 

The  cavitating  transfer  function  is,  however,  quite  different.  As 
our  touchstone,  the  slope  of  the  performance  curve  decreases  (in  magnitude) 
at  higher  frequency  and  the  inertance  part  practically  vanishes.  We  may 
mention  that  the  cavitation  number,  a  =  0.046,  results  in  a  condition  of 
extensive  cavitation  in  the  inducer  portion  of  the  pump.  There  is  attached 
blade  cavitation,  but  the  extensive  cavitation  seen  under  strobe  illumi¬ 
nation  is  due  primarily  to  bubbly-cavitating  tip-clearance  flows  (see 
Figure  41).  This  extent  of  cavitation  is  a  normal  condition  of  operation 
for  many  inducer  pumps.  We  see  some  change  in  the  customary  pressure 
rise-flow  characteristic  (z^)  due  to  cavitation.  But  there  are  more 
important  and  striking  changes  in  the  other  terms  of  the  transfer  function; 
the  term  z^,  for  example,  exhibits  a  large  change  with  frequency,  both  in 
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Figure  40  -  Test  Loop  for  Dynamic  Transfer  Function  Measurement 

(Ng  1976) 
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Figure  41  -  Transfer  Function  for  an  Axial  Inducer  Pump  at  a  Flow  Co¬ 
efficient  of  0.07  (Inlet  Value)  (The  value  of  the  total  pressure 
coefficient  is  about  0.50.  The  inducer  consists  of  four 
cambered  blades  having  a  tip  inlet  angle  of  7  degrees 
followed  by  12  more  heavily  loaded  tandem  blades. 

The  real  part  of  the  coefficients  are  shown 
solid  and  the  imaginary  dotted.  Two 
sets  of  data  are  shown,  fully 
wetted  and  for  a  cavitation 
index  o  =  0.046.) 
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the  real  part  and  in  the  imaginary  part.  The  conventional  cavitation  test 
at  this  operating  condition  shows  a  very  slight  reduction  in  output,  a 
finding  consistent  with  the  low  frequency  limit  shown  in  Figure  AO.  We  see 
under  dynamic  conditions  a  much  larger  change  in  output  pressure  with 
inlet  pressure  thai  occurs  in  steady  flow;  moreover,  there  is  a  pronounced 
phase  change  for  the  imaginary  part.  We  must  conclude  that  this  part  of 
the  cavitation  performance  is  complex  (i.e.,  having  real  and  imaginary 
parts)  and  is  a  strong  function  of  frequency.  Similar  deductions  flow 
from  the  terms  z2^,  the  compliance,  and  z^,  which  has  been  termed  by 
Brennen  as  the  mass-flow  gain  factor.  The  overall  impression  from  the 
experiments  summarized  in  Figure  41,  is  that  there  are  significant  dynamic 
effects  in  the  pump  transfer  matrix  and,  importantly,  that  these  effects 
cannot  be  estimated  from  quasi-steady  tests  or  attached  cavity  cascade 
theories  (see  Ng  and  Brennen  for  details). 

We  have  already  drawn  attention  to  the  bubbly  nature  of  this  particu¬ 
lar  flow  in  Figure  42.  With  this  in  mind,  Brennen  (1978)  has  formulated 
a  dynamic  model  of  the  flow  through  inducer  pump  blade  channels  based  upon 
a  bubbly  two-phase  flow  concept.  Although  this  is  empirical  in  some 
respects,  it  reproduces  most  of  the  effects  of  frequency  and  cavitation 
number,  seen  in  Figure  41,  in  a  remarkable  way.  Inevitably,  again,  the  void 
fraction  depends  upon  some  initial  (but  unknown)  microbubble  population. 

Transfer  function  representation  of  pump  dynamics  is  now  becoming 
more  common.  We  have  alluded  to  the  work  of  Fanelli,  and  Black  and  Santos 
who  do,  in  fact,  use  representations  essentially  the  same  as  Equation  (39) 
for  the  governing  and  control  of  hydraulic  turbines,  and  in  analyzing 
self-excited  oscillations  in  pump  systems,  respectively.  These  later 
"auto"-oscillations  (distinct  from  the  Pogo  problem)  may  lead  to  damaging 
levels  of  pressure  fluctuations.  The  onset  of  such  oscillations  in  a 
dynamically  characterized  hydraulic  test  loop  has  been  shown  recently  to 
depend  on  the  measured  pump  transfer  function  (Braisted  and  Brennen  1978) . 
More  importantly,  Brennen  (1978)  has  argued  that  the  transfer  function 
itself  cannot  be  represented  simply  by  a  collection  of  passive  elements 
such  as  compliance,  resistance,  and  inertance,  i.e.,  the  conventional 
lumped  parameters  often  used  in  analysis  of  mechanical  linear  systems. 
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tation  in  the  Pump  of  Figure  41  at  the  Cavitation  Index  of 


We  conclude  here,  our  discussion  of  cavitating  internal  flows;  the 
field  remains  extremely  active  because  of  the  thermo-hydraulic  problems 
associated  with  nuclear  power  plants.  We  have  stressed  application  of 
these  internal  flows  to  pumps  because  of  the  relative  simplicity  of  the 
representation  and  because  of  the  inherent  importance  of  the  subject.  As 
we  have  seen  there  are  somewhat  analogous  problems  in  naval  hydrodynamics. 
Some  of  these  are  due  to  geometric  confinement,  such  as  in  the  testing  of 
unsteady  cavitating  flows  in  water  tunnels,  and  others  have  their  common 
origin  in  the  liquid  environment  itself. 

SUMMARY 

It  has  been  nearly  ten  years  since  "Cavitation  State  of  Knowledge" 
has  appeared.  Since  then,  many  new  lines  of  research  have  appeared  and 
some  older  ones  refined.  "There  is  now  a  much  clearer  understanding,  than 
formerly,  of  the  role  of  the  viscous  boundary- layer  flow  on  smooth  bodies 
in  cavitation  inception.  Similarly,  the  necessity  of  carefully  describing 
the  form  and  extent  of  cavitation,  both  at  its  onset  and  developed  states 
for  the  formulation  of  realistic  physical  models,  seems  well  appreciated. 
New  experimental  techniques,  including  the  measurement  of  freestream  cavi¬ 
tation  "nuclei,"  have  appeared;  these  promise  to  put  laboratory  and  field 
studies  of  cavitation  on  a  new  quantitative  level.  The  presently 
available  cavitation- incept ion  scaling  concepts  are  considerably  improved 
and  more  physically  based  than  formerly.  With  all  these  new  advances  it  is 
still  not  possible  to  predict  with  certainty  the  location,  form,  and  onset 
conditions  required  for  cavitation  in  an  arbitrary  flow.  This  is 
particularly  so  for  unsteady  cavitating  flows.  These  flows,  in  addition, 
because  of  the  volume  change  of  the  cavitation  itself,  may  experience  large 
unsteady  interactions  with  neighboring  surfaces  or  flow  fields.  Initial 
attempts  at  sorting  out  these  mutual  effects  for  certain  well-defined 
hydraulic  systems  have  met  with  some  success.^  As  our  observational  and 
methodological  techniques  improve,  there  is  every  reason  to  believe  that 
our  understanding  of  unsteady  cavitation  phenomena  will  continue  to  im¬ 
prove,  just  as  it  has  for  cavitation- inception  processes. 
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